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In  this  dissertation,  the  output  feedback  stabilization 
problem  is  studied  from  a  more  general  framework.   For  any 

ex  <  B,  let  Snx    and  C^  p  denote  {(a^  a2 an)  t    Rn  |  xn  + 

Z j=la jx"  J    =    0  implies  Re (x)  <  a}  and  { (a± , a2 , . . . ,an)  e  Rn | 
x   +   j=iajxn  D  =  0  implies  |x  -  (a+B)/2|  <  (B-a)/2)  respec- 
tively.  In  chapter  one,  we  illustrate  that  under  the  general 
setting,  the  output  feedback  stabilization  problem  of  a 
system  with  single  control  can  be  converted  to  a  problem  of 
examining  whether  a  given  affine  set  (depending  on  the 
system)  intersects  SQ.   Main  results  of  this  dissertation  are 
also  mentioned.   In  chapter  two,  we  obtain  easy  necessary 
and  sufficient  conditions  for  a  hyperplane  H  to  intersect 
s5-t  Wc  also  show  that  H  n  s"  *    <{.  if  and  only  if  II  n  RsJJ  *    <f> 
where  Rs"  =  {  (ax  ,  a2 ,  .  .  . , aR)  c  Rn |  all  the  roots  of  xn  + 
Zj=lajx  J    =  0  are  real,  distinct  and  strictly  less  than  a}. 
The  results  are  based  on  the  following  property:   "The  convex 
hulls  of  RSQ  and  S0  are  both  equal  to  the  positive  orthant 


of  R  •"   J"  chapter  throe,  we  characterize  hyperplanes  which 
intersect  CQ  „  for  arbitrary  fixed  a  <  6.   We  also  show  that 
c°(Ca,8^  is  an  °Pcn  simplex  in  Rn  generated  by  some  special 
vertices.   In  chapter  four,  we  investigate  the  output  feedback 
stabilization  problem  of  a  special  delayed  system.   Chapter 
five  contains  a  sufficient  condition  for  a  straight  line  to 
intersect  SQ  and  a  computer  aided  algorithm  for  determining 
the  SQ-stabilizability  of  a  given  straight  line.   Some  dis- 
cussion of  the  directions  for  future  research  is  included. 


CHAPTER  I 
INTRODUCTION 


Consider  an  unstable  continuous  linear  time  invariant 
control  system: 

(C)   X  =  Fx  +  Gu,  y  =  Hx 

where  F  e  Rnxn,  G  e    RnXm,  H  e  RlXn.   One  of  the  problems  which 
have  been  receiving  a  great  deal  of  attention  in  control 
theory  is  the  Output  Feedback  Stabilization  problem.   It  is 
to  determine  the  existence  of  an  output  feedback   u   =  Ky  or 
even  to  find  such  K  e  R™  r  so  that  the  closed  loop  system 
becomes  stable.   In  other  words,  it  is  to  find  a  matrix 
K  e  RmXr  such  that  all  the  roots  of  the  following  equation 
have  negative  real  parts. 

det(xl  -  I  -  GKH)  =  0 

Wonham  [1]  had  proved  the  following  theorem. 

Theorem  1.1.1.   For  a  control  system  defined  as  in  (C),  let 
H  =  I.   Then  the  pair  (F,G)  is  controllable,  namely,  rank 
(G,  l'G,  .  .  .  ,  Fn~  G)  =  n  (see  Kalman  et   a  1  .  Ill)  if  and  only  if 
for  every  choice  of  (a,f a-, . . . , a  )  c    Rn ,  there  is  a  matrix 
K  e  Rr"Xr  such  that  det(xl  -  F  -  GKH)  -  x"  +  }.n_  ^a  .  xn~ J  .   I ] 
Note  that  in  this  theorem,  physically,  the  assumption  H  =  I 
implies  that  all  the  states  are  measurable.   This  is  unlikely 
to  most  practical  situations. 

1 


A  considerable  amount  of  research  in  the  Literature  has 
been  devoted  to  finding  a  controller  designed  from  the  avail- 
able state  measurements  to  stabilize  the  system.   One  approach 
is  to  obtain  estimates  of  the  inaccessible   states  from  a 
Luenberger  Observer  (Luenberger  [1],  [2],  [3])  and  use  these 
estimates  along  with  the  available  state  measurements  to 
construct  a  feedback  control.   Along  the  same  line,  Brasch 
and  Pearson  [1]  showed  that  arbitrary  pole  placement  can  be 
obtained  by  adding  a  compensator  to  the  system.   Both  are 
powerful  methods  for  stabilizing  a  system;  yet,  one  pays  the 
price  of  increasing  the  dimension  of  the  system.   Hence,  it  is 
still  desirable  to  know  whether  a  constant  gain  feedback 
can  be  designed  directly  from  the  available  state  measurements. 
Most  of  the  efforts  of  the  researchers  such  as  Miller,  Cochran 
and  Howze  [1],  Luus  [1],  McBrinn  and  Roy  [1],  and  Sirisena  and 
Choi  [1]  were  devoted  to  the  study  of  numerical  schemes  for 
finding  the  feedback  gain  matrix  k.   In  a  paper  by  Anderson, 
Bose  and  Jury  [11,  the  problem  of  output  feedback  stabilization 
was  related  to  the  decidability  problems  (see  Jacobson  [ID 
and  these  authors  showed  how  the  decision  methods  of  Tarski 
[1]  and  Seidenberg  [1]  can  be  applied  to  solve  the  stabilization 
problem. 

In  this  dissertation,  we  shall  look  at  the  output  feedback 
stabilization  problem  from  a  more  general  point  of  view 
described  in  the  following. 


Definition  ];   Given  a  set  K ,    an  n-di monsional  linear  output 

feedback  syr>tern  with  parameters  in  K    is  defined  to  ho  an 

ordered  pair  (tt,K)  ,  where  u  is  a  mapping  from  K    into  RnXn  x 

nxm     mxn  . 

R     *  R     for  some  fixed  m,  which  associates  with  each  para- 


meter K  e  K       a  triple  (A(K),  D(K),  C(K))  with  A(K)  c  r"*"1, 

For  brevity,  in  the  sequel  we 


(K)  c    RUXm    and  C(K)  e    RmXn , 


simply  say  "a  system  (it  ,  K)   with  parameters  in  K .  " 

For  a  fixed  m  and  n,  let  x  be  a  mapping  from  RnXn  x  Rr 

mxn  .      n   ,  .  , 
R     into  R   which  assigns  each  (A,B,C)  a  vector  (a, ,a9,.. 

an)  £  R   whore  a.'s  satisfy  x"  +  Zn=1a.xn_:]  =  det(xI-A-BC) 

Let 


Sn  =  {(a., a-,..., a  )  e    Rn |  xn  +  I 

a  1   2       n       ' 


.  .a.x 
D  =  l  1 


n-j 


=  0  implies 


Re(x)  <  a} 

RS"  =  { (a, ,a_, . . . ,  a  )  e  Rn |  all  the  roots  of  xn  +  Zn  na.xn~i 
u       iz       n  l=li 

are  real,  distinct  and  strictly  loss  than  a} 


Definition  2:   For  a  fixed  a  e  R,  an  n-d imensional  linear 

output  feedback  system  (t\,K)       with  parameters  in  K,     is  said 

to  be  S^-stabilizable  if  xT,(^')  n  Sn  *  <J, .   A  system  (n,k')n  with 

parameters  in  K  is  said  to  be  RSn-stabilizable  if  xtt(K)  n  RSn  *  c. 

a  a 

If  one  lets  \  =  R     and  defines  n  (K)  =  (F,G,KH)  for  each 
K  e  K   where  F,G,H  are  the  matrices  given  in  (C),  then  xTT(^) 
is  equal  to  the  coefficients  of  the  characteristic  polynomial 
of  (C)  with  output  feedback  u  =  Ky .   Note  that  (C)  is  stabiliz- 
able  if  and  only  if  x11^')  n  sn  *    '!' /  or  equiva  lent  ly ,  by 
definition,  (tt  ,K)       is  S„-stabi  li  zable .   Suppose  H  =  I  and 
(F,G)  forms  a  completely  controllable  pair.   Under  the  same  it 


as  above,  we  observe  that  the  wel]  known  result  (theorem  1.1.1) 
of  Wonham  [1]  means  then  >; n  ( '0  =  Rn . 

For  any  system  (C)  with  single  control  i.e.,  m  =  1,  it 

was  shown  that  (see  Byrnes  [1],  Karnen  ilj) 


det(xI-F-GKH)  =  det(xI-F)  +  kH  adj(xI-F)G 


lxr 


( i .  i ; 


With  K  =  R     and  the  same  it  as  above,  the  relation  (1.1) 

implies  that  vtt(^)  =  x^(0)  +  KM  where  M  is  a  r  x  n  matrix 

depending  on  F,G,H.   Hence,  x"n  (K)     forms  an  affine  set  in  Rn . 

Now  consider  the  class  of  all  systems  (7i,K)n  for  which 

XTt(K)  forms  an  affine  set.   This  is  a  large  class  which  includes 

all  systems  (C)  with  single  control  as  the  discussion  in  the 

above  shows.   Observe  that  for  a  fixed  a  e  R  if  one  can 

characterize  affine  sets  which  intersect  Sn,  then  Sn-stabili- 

oc         a 

zability  of  a  system  (v,K)n    for  which  v (K)     forms  an  affine  set 
is  characterized.   In  particular,  the  problem  of  output  feed- 
back stabilization  for  systems  with  single  control  is  resolved. 
As  a  first  step  in  this  direction  we  consider  mainly,  the 
case  of  hyperplanos  in  this  dissertation. 

Chapter  2  is  devoted  to  the  study  of  characterizing  hyper- 
planes  which  intersect  S   and  RSn.   We  prove  that  a  hyperplane 


n  +  I.    .A.  a. 

0     ]=1  j    j 


0  }  i  n  t  e  r  s  c  c  t  s  S  _.  i  f 


H  =  {  ( a  .  ,  a  .  ,  .  .  .  ,  a  )  e  Rn  I  A 
12       n       ' 

and  only  if  !i  n  R      where  R   is  the  positive  orthant  of  Rn,  i.e., 
{ (a1,a2, . . . , an)  t  R  |  a .  >  0}.   Furthermore,  it  is  shown  that 
H  n  Sn  *  <|i  is  equivalent  to  II  n  Rs"  ^  $.       We  obtain  also  the 


criterion  which  states    H  n   SA  *  $  if  and  only  if  A. A.  <  0 

0  !   j 


for  some  i  /  j,  i,j  =  0  ,  1 ,  2  ,  .  .  .  ,  n  .   The  simplicity  of  those 
properties  results  from  the  following  lemma: 


"The  convex  hulls  of  sJJ  and  Rsj]  are  equal  to  the  positive 
orthant  r"  of  Rn  i.e.,  co(Sq)  =  co(RsJj)  =  r"." 


The  author  would  like  to  acknowledge  that  for  the  case 
of  n  =  4,  the  lemma  was  first  shown  by  Professor  Popov.   Based 
on  the  idea  from  that  proof,  this  author  extended  the  lemma 
for  arbitrary  n.   Also,  the  proof  for  general  n,  was  greatly 
simplified  by  a  suggestion  from  Professor  Popov. 

The  results  mentioned  above  are  generalized  to  characterize 

all  hyperplanes  which  intersect  Sn  or  RSn  for  an  arbi trarv 

a       a  J 

fixed  a  e    R.   Observe  that  for  any  system  (;r,K)n  with  param- 
eters  in  K,    there  exists  always  a  -:  R  such  that  ytt(K)  n  Sn  *  & 

a 

From  the  point  of  view  of  stability,  it  is  desirable  to  have 

that  a  "as  negative  as  possible"  to  reduce  the  transient  time. 

However,  if  >;tt(/C)  forms  only  a  hyperplane  in  Rn ,  then  one  can 

not  expect  the  relation  X*  (K)     n  Sn  *    ,|  to  hold  for  any  a  e  R. 

We  show  that  given  a  hyperplane  H  the  value  inf{a|  H  n  Sn  *  4} 

a 

is  equal  to  the  negative  of  the  maximal  real  root  of  a  specially 
constructed  polynomial.   Some  computation  aspects  such  as  how 


,a  )  c  H  n  S  .  if  H  n  S11  *    <],,  are  dis- 
n  a  a 


to  find  a  (a, , a~ , 

cussed.   An  example  chosen  from  Anderson  et   al.  [11  is 

included . 

Given  a  discrete  linear  time  invariant  control  system 


(5 


\  +  l 


FxR  +  Gu,  yv 


k  ~  Hxk 


(S)  is  stnbi ] izablc  with  output  feedback  if  and  only  if  there 

exist  K  c  R     such  that  the  modulus  of  all  the  roots  of 

det(xl  -  F  -  GKH)  =  0  are  less  than  one.   Given  a    <    B,  let 

Cr,  ft  =  Ua,  ,a7,  .  .  .  ,  a  )  e  R  |  x   +  In  ,  a.xn_:]  =  0  implies 
ex ,  p       Ll               n       '        1 =  1  1 

|x  -  (6+a)/2  |  <  (6-a)/2) 

Re"  p  =  {(alfa0f...,a  )  e  Rn |  all  the  roots  of  xn  +  Zn  na.xn~J=0 
a ,  p       VI  n  j  =  l  j 

are  real,  distinct  and  strictly  between  a  and  (3 }  . 


By  taking  K    =  R     and  tt(K)  -  (F,  G,  KH)  for  each  K  £  K,  one 

can  regard  (S)  with  output  feedback  as  a  system  (tt,K)   with 

parameters  in  K.   (S)  is  stabilizable  by  an  output  feedback 

n 


if  and  only  if  xn(^)  n  c 


1,  1 


The  main  objective  of 


chapter  3  is  to  find  a  necessary  and  sufficient  condition  for 

a  hyperplane  to  have  nonempty  intersection  with  C_ ,  ,.   For 

the  reason  of  mathematical  interest,  we  obtain  necessary  and 

sufficient  conditions  for  a  hyperplane  to  intersect  C   o  for 

general  a    <    $.       It  is  interesting  to  point  out  that  given 

a  <  B,  the  C^  c  and  RC  0    are  both  equal  to  the  open  n  simplex 
a. ,  p        a ,  p  -1  '         : 

in  R   aenerated  by  the  n+1  vertices  (q ., ,q ._,... ,q .  )  c  R 

l  =  0  ,  1 ,  .  .  .  ,  n  where  x   +  I  .  _•,  q  .  .x     =  (x-cx)     (x-p)  . 

In  chapter  4,  we  consider  the  stabilization  of  a  special 
kind  of  system  with  delay: 

(P)   x(t)  =  ax(t)  +  bx(t-l)  +  u,  y  =  cx(t)  +  dx(t-l),  u  =  ky 


With  a,b,c,d  fixed,  one  wishes  to  find  k  e  R  such  that  all  the 
roots  of  X  +  (a+kc)  +  (b+kd) e  '  =  0  have  negative  real  parts. 
Given  a  e    R,  let 


V 


{(a,b)  t     Rn|  X+n+bc  ?  =  0  implies  RcX  <  a} 


This  stabilization  problem  turns  out  to  be  a  problem  of 
determining  whether  a  given  straight  line  i nterseets  PQ.   We 
obtain  a  necessary  and  sufficient  condition  for  a  straight 
line  to  intersect  V       for  general  a  c    R.   It  is  interesting  to 
note  that  for  any  given  hyperplane  II,  the  number  inf{a|  II  n  Sa 

*  <j>}  is  always  bounded  from  below  (theorem  2.3.1).   However, 

2  2 

in  this  delay  case,  there  are  lines  H  c  R   (hyperplane  in  R  ) 

for  which  the  number  inf{cx|  H  n  PQ  *  <M  =  -«  (theorem  4.2.4). 

One  of  the  problems  concerning  the  stabilization  of 
systems  with  delay  is  that  "Can  one  stabilize  a  system  with 
delay  by  a  feedback  without  delay?"   In  the  end  of  the 
chapter  4  we  show  with  an  example,  that  this  can  not  be  done 
in  general,  despite  the  fact  that  this  can  be  done  for  (V) . 

Despite  all  the  easy  criteria  for  determining  different 

kinds  of  stabilizabilities  obtained  for  hyperplanes  in  chapters 

2,  3  and  4,  we  are  unable  to  find  an  easy  necessary  and 

sufficient  condition  for  affine  sets   with   dimension  strictly 

lower  than  n-1  to  intersect  s".   As  an  attempt, in  chapter  5, 

we  obtain  a  sufficient  condition  for  a  straight  line  to 

intersect  Sn.   Also,  a  computer  aided  algorithm  is  proposed 
a 

for  determining  the  Sa~s tabi lizabi li ty  of  a  straight  line  . 
Some  directions  for  future  research  are  discussed. 


n  CHAPTER  II 

Sa-STABILIZABLE  AND  RS"-STABILIZABLE  HYPERPLANES 

For  a  fixed  n,  let  (tt,K)   be  a  n  dimensional  system  with 

parameters  in  K.   Recall  that  K  is  a  set  of  parameters  and  tt 

is  a  mapping  from  K  into  RnXn  *  RnXm  x  RmXn  for  some  fixed  m, 

which  associates  each  parameter  K  e    K,    a  triple  denoted  by 

(A(K),  B(K),  C(K)).   Let  x  be  the  mapping  from  Rnxn  x  Rn*m  x 

RmXn  into  Rn  which  assigns  each  (A,B,C)  e  Rn*n  x  RnXm  x  Rmxn 

a  vector  (a..  ,  a0  ,  .  .  .  ,  a  )  e  Rn  where  xn  +  Zn  ,a.xn"]  =  det(xl  - 
1      Z  n  ]=1  j 

A  -  BC) . 

The  main  purpose  of  this  chapter  is  to  characterize  the 

S  -stabilizability  and  RS  -stabilizability   of  systems  (n,K)n 
a  a  J 

for  which  xv  {k)     forms  a  hyperplane  in  R  .   In  other  words,  one 

wishes  to  characterize  hyperplanes  which  intersect  Sn  and  RSn. 

a        a 

Recall  that   for  a  given  a  e  R, 


n 

Sa  =  { (ai'a2' • ' • 'an}  e  rR!  x"  +  l    a-xn_D  =  °  implies  Re(x)<a] 

j  =  l  D 


RS"  =  {  (a, ,a„, . . . ,a  )  e    Rn I  all  the  roots  of  xn  +  I    a . xn  ^  =  0 

1   2  j=i  3 

are  real,  distinct  and  strictly  less  than  a). 

Definition :   Given  a  e  R,  a  hyperplane  H  <  RU   is  said  to  be  Sn- 

stabilizable  if  H  n  S  *    § .  A  hyperplane  H  is  said  to  be 

RS^-stabilizable  if  H  n  RSn  *  <f> . 

a  a 


The  technique  which  we  use  in  here  is  based  on  a  simple 
property  that  a  hyperplane  in  R   intersects  an  open  connected 
set  in  R   if  and  only  if  the  hyperplane  intersects  the  convex 
hull  of  that  open  connected  set.   In  lemma  2.1.5,  we  prove 
that  co(RSq)  =  co(S  )  =  R+.   By  an  invertible  affine  transfor- 
mation T    and  lemma  2.1.5,  we  obtain  the  descriptions  of 

co(S  )  and  co(RS  ) .   Section  2  contains  some  criteria  for 
a  a 

determining  the  S  -stabilizability  of  a  hyperplane.   We  show 

that  the  S  -stabilizability  and  RS  -stabilizability  are  eauiv- 

alent  for  hyperplanes.   In  particular,  we  find  that  a  hyperplane 

is  S„-stabi lizable  if  and  only  if  the  hyperplane  intersects 

the  positive  orthant.   Finally,  in  section  3,  the  value 

inf{a|H  n  S   *  <j>}  for  a  given  hyperplane  H  c  R   is  computed 

and  some  consideration  of  actually  finding  a  point  (a. , a~ , . . . , 

a)eHnS       *    i>    is    included, 
n  a         r 

2.1   Convex  Hulls  of  Sn  and  RSn. 

— . , . _ _a q 

Let  us  beain  with  some  basic  properties  of  S   and  RS 

^  f   f  a        a 

which  will  be  needed  later.   Recall  that  S   =  {£.  ,  a.e.  e  R  I 

a      i=l  i  l      ' 

x  +  £.  na.x  J  implies  Re(x)  <  a}  and  RS  =  {£  .  ,a-e-  c  R  I 
all  the  roots  of  x  of  x  +  E._,ax  -1  =  0  are  real,  distinct 
and  x  <  a } . 

Lemma  2.1.1.   (i)   RS   and  S   are  open  connected  sets  in  R 
a       a 


for  any  a  c    R; 

(ii)   Sn  c  s"  if  and  only  if  a    <    6 ; 

u  p 

(iii)   RSn  c  RS"  if  and  only  if  ex  <  6; 

ex  p 

(iv)   RSn  c  Sn  c  r"  for  every  a  <  0. 

a  a     + 


The  proof  of  ( i )  will  be  presented  in  appendix   I 
the  others  follow  straightforwardly  from  the  definitions  of 
Sq  and  RS^  together  with  the -fact  that  Sn  c  Rn. 

Since  Hermite  [1],  enormous  research  had  been  devoted  to 
the  study  of  the  set  s£.   Algebraically,  SnQ    can  be  characterized 
by  the  well  known  criteria  of  Routh  and  Hurwitz  or  Lidnard 
and  Chipart  (Gantmacher  [1]),  in  terms  of  system  of  nonlinear 
inequalities.   From  these  criteria,  it  follows  that  Sn  is 
subset  of  R+.   Since  the  roots  of  a  polynomial  depend  contin- 


uously on  the  coefficients,  the  boundary  of  s"  consists  of 
coefficients  of  polynomials  which  have  either  a  zero  root  or 
a  pair  of  purely  imaginary  roots.   Indeed,  the  D-decomposition 
method  of  Yu  Naimark  says  that  the  hypersur faces  H,  =  {(a  ,a 

•••'an)  e  ^1  an  =  0}  and  H2  =  {(al'a2 an)  e  R"  I  x"  + 

j=lajx     "  °  has  a  pair  of  purely  imaginary  roots}  divide 
the  coefficient  space  Rn  into  n+1  disjoint  regions;  and  the 
polynomials  with  coefficients  chosen  from  the  same  region  will 
have  the  same  number  of  roots  with  positive  real  parts.   The 
example  below  illustrates  the  case  n  =  2. 


Example:   For  n  =  1,  SJ  =  {(a-^a^l  a± 


0,  a2  >  0}  and 


RS0  ~  ((a1,a2)|  a1  >  0,  a2  >  0,  a1    >    4a^}.   On  the  curve 

2 
al  "  '4a2'  thc  Polynomials  will  have  identical  roots.   In 

figure  2.1,  p  represents  the  number  of  roots  with  positive 
real  parts. 

The  sets  SQ,  for  a  (    R,  as  exemplified  by  s"   are  deter- 
mined by  collections  of  nonlinear  inequalities  for  a.'s.   In 


^w< 


Figure    2 . 1 


general,  for  a  given  hyperplane  II  <~    R  ,  it  is  not  easy  to 

determine  if  H  is  S  -stabilizablc  by  considering  directly  the 

intersection  of  H  with  S  .   However,  because  S   is  open 

a  a 

connected  in  R  ,  the  following  lemma  shows  that  one  can 

consider  instead  the  intersection  of  H  with  co(Sn),  a  set 

a 

which  is  larger  and  more  regular  than  S1  itself. 

Lemma  2.1.2.   Let  H be  a hyperplane  in  Rn ,  and  let  S  c  Rn  be  an 
open  connected  set.   Then  H  n  S  *  (J)  if  and  only  if  H  n  co(S)  *  it- 
Proof:   (=>)  .   Trivial,  since  S  c  co(S). 

(**=)  .   Let  E,  and  E_  be  the  open  half  spaces  corresponding 
to  the  hyperplane  H.   Assume  the  contrary  that  H  n  S  =  <J>. 
Then,  by  the  hypothesis  that  S  is  open  and  connected,  we  have 
either  S  c  E ,  or  S  c  E  .   Without  loss  of  generality,  let  us 
assume  S  £  E,  .   Since  E,  is  convex  and  S  c  E, ,  it  follows  that 
H  n  co(S)  c  H  n  co(Ej  =  H  n  E,  =  (j> .   This  contradicts  the 
hypothesis  H  n  co(S)  *  (J).   Therefore,  one  concludes  that 
H  n  S  *  <J>.  0 

Using  this  lemma,  the  S  -stabi lizable  hyperplanes  can  be 

characterized  as  those  which  intersect  co(S  ) .   However,  if 

a 

one  wishes  to  take  the  full  advantage  of  this  characterization 

of  S  -stabilizablc  hyperplanes,  one  must  know  what  co(S  ) 
a  j  i   r  a 

looks  like.   This  will  be  the  aim  of  the  rest  of  this  section. 
First,  we  shall  prove  the  main  lemma  mentioned  in  the  intro- 
duction (i.e.,  co(S  )  =  co(RS  )  =  R  ) .   Then,  we  shall  use 

this  lemma  to  obtain  descriptions  of  co(S  )  and  co(RS  )  for 

a  a 

arbitrary  a  e  R. 


1  i 


In  order  to  prove  the  main  lemma,  since  co(Sf),  co(RSn! 
and  R,  are  open  convex    sets,  it  suffices  to  show  that  their 
closures  are  equal  i.e.,  co(Sq)  =  co(RSq)  =  R+  (see  e.g., 
Rockafellar  Ll],  p.  45) .   This  is  what  we  establish  in  the 
following  lemmas. 


Lemma  2.1.3.   The  origin  6  =  (0,0,...,0, 


in  co(RS„) 


,  Oe  .  is  contained 

1  =  1   1 


Proof:   Let  6  be  a  positive  number.   Consider  the  following 
polynomial 


n  (x+j6)  =  x  +  i   p. 

j=l  j=l  D 


5)x 


n-i 


Observe  that  for  each  j  ,  P.  (6)  =  c..o-'  where  c.  is  a  positive 

real  number.   Therefore,  for  any  e  >  0,  one  can  find  6  >  0 

-v      -1/2 
such  that  P-(o)  <  n  '     e  for  every  j  =  l,2,...,n.   Note  that 

n 


II  ._  ,(x  +  j6)  =  0  has  only  stable  roots  x  =  -j<5,  j  =  1,2,..., n 
Thus,  by  the  definition  of  RSQ, 


Z    P  .  (  6  )  e  .  e  RS  _ 
=  1  3     1      0 


(2.1) 


Also,  one  has 


P  .  (6)e  •  -  I    Oe. 


1  =  0  3 


j  =  l 


(  Z    P  •  (  6  )  ) 
i=l  D 


1/2 


(  I 
j  =  l  n 


2,  1/2 
c  )     =  c 


By  (2.1),  it  fellows  that  Y.  .     .P. 
1  3  =  1  ] 


P  .  ( 6  )  e 


Hence,  0 


co  (  RS 


N  (0)  n  RS„. 
c         0 


Let  {e.}.  ,  be  the  standard  basis  of  R 
1  1  =  1 


1  '1 


Lemma  2.1.4.   For  any  r  >  0  and  any  positive  integer  k,  the 
point  re,  is  contained  in  co(Rs") . 

Proof:   In  this  proof  j  will  be  used  to  represent  positive 
integers.   Take  r  and  k  as  specified  in  the  lemma. 

For  arbitrary  q  >  0  and  6-  >  0,  i  =  1,2, ...,n  where 
<5j_  *  6.  if  i  *  j,  let  P  .  (q,  6.  ,  62 ,  -  -  -  ,  &    )  be  the  polynomial  of 
q  and  fi.'s  (also  written  P.(q,6-))  obtained  from  the  following 
polynomial  expansion 


k 

II  (x+q  +  6 
i=l 


n   (r+6.)  =  xn  +  Z    P  (q,6   5     .,6  . 
i=k+l  j=l  D  n 


rn'3 


Note  that  for  every  j  <  k ,  the  highest  degree  of  variable  q 
in  P.(q,5.)  is  j.   In  particular, 


;(q,<V  =  qK  +  Pk(g,6.; 


2.2) 


The  degree  of  q  of  P,(q,6.)  is  at  most  k-1.   Fix  5.  =6.  >  0, 
i  =  1,2,. ..,n.   For  every  e  >  0,  choose  q„  >  0  large  enough 
such  that  for  all  j,  1  <  j  <  k-1 


T    P  .  (qn  , 6,  ,6- ,  ....  6  )  <  i— 
k   j   '  0  '  1 '  2  '    '  n     v  n 


2.3! 


and 


"T  Vq0'*l'?2 V 


i'n 


(2.4) 


as  wc ] 1  as 


max{ 2 , r }  < 


^0 


(2.5) 


This  is  possible  because  the  degrees  of  a  in  P  (n,6  )  are  loss 

J     i 

than  or  equal  to  k-1  for  all  j  •  k-1  and  the  same  property  is 

true  for  P.  (q,6.).   Since  for  any  j  >  k,  P  .  (qn  ,  0  ,  0  ,  .  .  .  ,  0)  =  0 

and  since  P . (q_ , 6  , 6 _,..., 6  )  are  polynomials  of  6.'s  with 
3^012       n       l       J  i 

positive  coefficients,  one  can  find  0  <  6 .  <  6 . .  i  =  l,2,...,n, 
t  l     l  '       '  '    '  ' 

where  6.  *    6.  if  i  *  j.  such  that  for  each  i  >  k 


He  Pj(^0'6l'62 V  <    /n 


(2.6) 


Also,  since  for  each  i,  6-  <  6-,  by  the  property  that  P.(q,5- 
j  =  1,2,... ,n  have  positive  coefficients,  one  obtains  for 
each  j  <  k-1 


"I  Vq0'6l'62' 


;  )  <  -f  P.  (q_,  6.  ,6_,  ...  ,6  )  <  r- 

n      k   •]  J0   1'  2       n     /n 

cln 


2.  7! 


and  similarly, 

tVV{1{2 V     <    "I  ?k(q0'*l'*2 V     </5 

qo  qo 

Now,  consider  the  following  polynomial  equation 
k  n  n        „ 


(2.8) 


I; 


(x+q0  +  6i)  .  J   (x+6^  =  xn  +   E  P  (q(),51,62 6n)xn^ 

3  =  1  3=k+l     "  3  =  1 


It  has  only  stable  roots  x  =  - (q  +6 • )  <  0,  i  =  l,2,...,k  and 

n 
x  =  -6.  <  0,  i  =  k+l,...,n.   This,  by  the  definition  of  PS„, 

implies  that  the  point 


E  P  (q0,61,62 6n)e\ 

-,  =  i  J 


2.9: 


belongs  to  KSQ.   Hence,  using  Lhe  relation  RS"  r    eo  (  us'.')  ,  we 
obtain 


z   e  co(RS') 
e         0 


.  n 


From  Lemma  2.1.4,  6  =  ^^Oe.  .:  co(Rs").   Since  co(Rs")  is 


convex,  the  set  co ( RS  )  also  contains  the  line  segment  with 


k, 


the  origin  and  z£  as  end  points.   In  particular,  (r/q*)z   e 
co(RsJJ),  since  0  <  r  <  qjj .   By  (2.9)  and  (2.2),  we  have 


k    k  *e 


k-l 


=  {  Z  ("T  P^n'M   +  (-T7  P(q„,6,)) 


j=l  q 


k   j  ^0'  i 


I        (4-  P.  (q   6  ))2}1/2  (2.10) 

k  +  1  a.   -1   u   x 


j=k+l  q 


k  ->• 
We  now  prove  that  ||  reR  -  (r/qQ)  z  ||  <  c.   Let  R1IS  denote  the 

right  hand  side  of  (2.10).   Then,  by  (2.6),  (2.7),  and  (2.3) 

one  has 

2     2  2 

RHS  <  {(k-l)-—  +  —  +  (n-k)—  }1/2  =  c 
n     n         n 

Therefore,  we  have  shown  that  for  every  e  >  0,  there  exists 
a  point  -£  z£  e  N£(re}.)  n  co(RSq).   Hence,  re^  e  co(RSg).G 


q 


Now  we  state  the  main  lemma 


Lemma  2.1.5.   Tor  every  positive  integer  n,  co ( RS  )  =  co(s") 


Proof:   First,  observe  that  R™  ={E"   r,  (re,  )|r  >  0, 

+       K— 1  K     K 

Zk=lrk  =  lr    rk  >  °*'   since'  by  lemma  2.1.4,  re,  e    co(Rs") 


1  / 


for  every  r  >  0  and  every  k  =  ],2,...,n.   Hence,  it  follows 
that  r"  «:  co(co(RSq))  =  co(RSJJ) 


Conversely,  by  lemma  2.1.1 

lis       r'h      ml        r  Zn 


we  have  co(RSQ)  c  co(SQ)  c  co(R+)  =  R+ .   Therefore,  R/_  = 


co(RSn)  =  co(S,J.   Since  for  any  n  dimensional  convex  set 
S  c  Rn,  interior  (S)  =  interior  (S)  (see  Rockafellar  [1], 
p.  46)  also,  since  co(RS„),  co(S„)  and  R   are  open,  one  con- 
cludes that  co(RSq)  =  co(sJJ)  =  r".  G 

Before  we  set  out  to  find  RSn  and  Sn,  let  us  mention  the 

a       a 

following  remark. 

Remark:   If  T  is  an  invertible  affine  transformation  of  R  , 
then  for  any  S  £  R  , 


T(co(S)  )  =  co(T(S)  ) 

This  follows  easily  from  the  following  facts: 

(1)  Every  point  z  e  co(S)  can  be  represented  as  a  con- 
vex combination  of  points  from  S,  i.e.,  z  =  Z.k.s.  with  s.  e  S 

and  k.  >  0,  E.k.  =  1; 
l        11 

(2)  T(E.k.s.)  =  E.k.T(s-),  since  T  is  invertible  affine. 
K  ill      11    l 

Using  this  remark,  one  observes  that  if  there  exists  an 

invertible  affine  transformation  T   on  R   such  that  T  (S  )  =  S 

o.  ex   u 

then 


T  (co(Su)  )  =  co(T  (S")  )  =  co(S") 

a      a  a      a  U 


Hence,  the  convex  hull  of  Sn,  namely,  co(S  )  will  simply  be 

the  set  of  all  points  (a,,a2,...,a  )  c    Rn  whose  image  under 

T   belonqs  to  R , . 
a       -       + 


Such  a  transformation  does  indeed  exist,  as  we  will  show 
in  the  next  lemma. 

For  any  fixed  a  c    R,  let  Tq  be  the  transformation  on  Rn 
such  that  for  every  ( a,  , a_ ,  .  .  . , a  )  e  Rn , 

Vai'a2 an)  =  (a\'*2 ln] 

where  ai ' s  are  obtained  from  the  following  polynomial  expansion. 

U+a)n  +  I    a.(x+a)n_j  =  x"  +  I    a.xn_j  (2.11) 

j=l  J  j=l  ^ 

We  remark  that  this  transformation  was  considered  in  a 
paper  by  Hughes  and  Nguyen  [1],  and   the  explicit  expression 
of  a.  was  also  given  there.   For  completion  we  include  the 
proof. 


Lemma  2.1.6.   Tf   is  an  invertible  affine  transformation  on  Rn 
such  that 


T  (RSn)  =  RS"   and   T  (Sn)  =  s" 
a    a       0         a   a      0 

Proof:   (i)   T   is  invertible  and  affine. 
a 

We  first  observe  that,  by  Taylor's  theorem,  every  polyno- 
mial p(y)  =  y   +  E._  -,  a  .  y   -1  can  be  expressed  as 

n   o(j)(a)      i 
p(y)  =  I      ^-^--i^Cy-a)3  (2.12) 

j  =  0   J- 

Thus,  by  substituting  y  =  x  f  a,  (2.12)  becomes 

n         „  n    (n-i )  ,  . 

I    a.(x+a)n_:J  =  P(x+a)  =  T,      P.  \     a)    xn^ 
j  =  0  3  D=0   (n-D)l 


\'J 


with  a   =  1 
1  s  j  s  n 


] 


Tn  view  of  (2.11),  it  follows  that,  for  each  j, 

(n-j) 


P'   ""  (a) 


■ "  - 1  x  _  n  - 1- 


i  =  0  J 


(2.13) 


To  show  that  Tq  is  affine,  let  U|,a' ,a^)  and  (a" 

a2'---'a^  bG  arbitrary  two  points  of  Rn .   Then,  for  each  i, 

1  <    i  <  n  and  for  0  <  A  <  1, 

a.  =   S  (""Jlo^^Aa!  +  (l-X)a'.']  -  X[  E  (n~f)  a  j_ia !  ]  + 
:    i=0  J~x                     x  x       i=o  J-1       a 


+  (1-A) 


(  .  .  )  a J   a  '.' 


i=0 


J-i 


Hence,  Tq  is  affine.   Represented  in  matrix  form, 


Val'a2' 


t  a 


(  a.,  -i  -,,  ..  .,an)QQ  +  ((J)af  (")a2, 


■1r«2. 


,n.  n> 
,  (  )  a  ) 


where  Q   =  (a 


n~i^]-i 


].  ■)  is  an  n  x  n  matrix  with  entries  q.  .  =  (.  7)a 
1 1  d  l  ]    v  ]  - 1 ' 

i  =  1,2,... ,n,  j  =  1,2, ...,n.   Observe  that  Q   is  an  upper 
triangular  matrix  v/ith  all  diagonal  entries  equal  to  one. 
Therefore,  Q   is  invertible.   This  implies  that  T   is 
invertible . 

(ii)   We  now  prove  that  T  (s|J)  =  sjj  and  T  (RSn)  =  Rs" 
From  the  way  in  which  Ta  was   defined,  it  follows  that  for 


every  number  >•   such  that  xn  +  Y  .     ,a.x"  ^ 
u  0     j  -■  1  j  0 

y0  =  xQ  +  a  satisfies  p(yQ)  =  yJJ  +  T^^-Y^3 


0 ,  the  number 
0  and  v  i  ce 


versa.   Notice  that  Re(xQ)  <  0  if  and  only  if  Re(yfl)  <  a. 
Hence,  (a^,a2 ,  •  .  .  ,  an)  e  Sa  is  equivalent  to  T  (a-,,  a2  ,  .  .  .  ,  a  )  = 


(ax,a2, 


-a, 


Sn.   In  other  words,  T  (S  )  =  S_.   The  saint 


proof  holds  for  T  ( RS  )  =  RS ' 
a    a       0 


2  0 


Remark:   For  every  l|?  ^.c.,  lft  q(x)  -  x"  ^  ;•:"   a.x    .   An 
elementary  calculation  shows  that  the  inverse  of  TQ  is 
defined  by 

^1(a1,a2 an)  =  (a^a^ ^n)Qa    4-  (  (J)  (-a)  ,  (J)  (-ex)  2  ,  .  .  . 

(£)(-«)") 

~      ~  ■     •  i  ~      .  ia  —  i .  .   .  i  —  i 

where  0   =  (q .  • )  is  an  n  x  n  matrix  with  q .  ■  =  (■  ■)   (-a) 

Furthermore,  T~   satisfies  the  relations,  TQ  (SQ)  =  Sq   and 
a     0       a 

Theorem  2.1.7.   Given  any  a  e  R,  the  convex  hulls  of  Sq   and 

RSn  are  both  equal  to  the  polygonal  convex  set  formed  by  the 
u 

intersection  of  the  following  half  spaces: 

l    [  (n_1)a:'~1  la.  >  0      i  =  1,2, n 

i=0  3-i 


Proof-   Let  T   be  the  transformation  defined  in  lemma  2.1.6. 
a 

The  theorem  follows  readily  from  the  discussion  preceding 
lemma  2.1.6  and  (2.13).   The  same  argument  holds  for  RSa-  D 

In  terms  of  the  usual  standard  basis  of  R  ,  one  may  write 

R+  =  {Ik=lrA  I  rk  >  0K   SinCG  Ta1(S0)  =  Sa  and  CO(S0)  =  R+' 
it  follows  that 

S1"=l("-!:l("°|i"krk'Vrk  "  °'  '2'14) 


r.n 


whore  for  each  i,  V.     .q,  .o.  is  the  ith  row  of  matrix  Q  . 
Hence,  we  obtain  another  way  of  expressing  co(S  )  . 


Theorem  2.1.8.   co(Sn)  and  co(RSn)  are  both  equal  to  the  open 

■ ■        a  a 


convex  cone  with  vertex  E  ._,(.)  (-a) -'e  .  ,  generated  by  the  r 
vectors  of  matrix  Q   =  (q^^);  explicityly,  q..  =  (._  . 


(  •  ■  )  ( -  a ) 


ow 


Note  that  the  entries  of  the  row  vector  ((-,)(-a),  ( ~)  (-ex)  ,  .  .  .  . 
(n)  (-a)n)  represent  the  coefficients  of  the  polynomial  (x-a) 
which  has  n  identical  roots  x  =  a. 


Example:   For  n  =  2  and  a  fixed  real  number  a,  co(S  )  =  co(RS 
1 a 

2 

{ (a, ,  a~)   2a  +  a,  >  0  and  a   +a,a+a_>0.  Itis  the 


ll'"2'  '  ""  '   1 
shaded  area  in  figure  2.1 


As  an  easy  generalization  of  theorem  2.1.8,  one  has 


Theorem  2.1.9.    Given  a  e  R,  suppose  T    is  a  set  of  complex 
numbers  which  satisfies 


{ x  I  x  <  a }  c  f  £  { z  I  Re  ( z )  <  u 


(2.15) 


Let  Il(r)  denote  {Zn=,a.e\|  xn  +  E^_  a.xn  :  =  0  implies  x  e 


r} 


Then,  co(IUD)  =  co(s")  =  co(RS"). 


Proof:   (2.15)  implies  that  co(RSn)  c  co(lKD)  £  co(S": 


By 


the  relation  co ( RS 


co(RS  )  =  co(S  ) . 


co(S  ),  one  concludes  that  co(iUI') 

a 


2.2      Criteria    of    S    -Stabi lizability    of    ilyperplanes 

a 

First,  let  us  prove  the  following  fact. 


Lemma  2.2.1.   Given  a  hypcrplnne  H  =  {  ( a  ,  a  , ,  .  .  .  ,  a  )  c  Rn  | 
A0  +  ^i  =  l^iai  ~  °^  tlie  following  properties  arc  equivalent 


(i)   H  n  R 


'!< 


(ii)   At  least  two  of  the  numbers  A.,  i    =    0,1.2, 
have  strictly  opposite  signs. 


,  n 


Proof:   (i)  =  (ii).   This  follows  easily  from  the  observation 

that  the  sum  A„  +  £._,A.a.  will  never  equal  zero  for  any 

(alfa2 , . . . ,an)  c  R+  if  all  A.,  j  =  0,1,2, ...,n  are  of  the 

same  sign. 

(ii)   *  (i).   Suppose  that  there  exist  i  *  j  such  that 

A^A.  <  0.   Without  loss  of  generality,  let  us  assume  A„  >  0 . 

Let  A  =  {j|  A.  <  0}  and  A  =  { j I  A.  >  0}.   Then,  H  n  Rn  *  A  , 
3  3  + 

since  H  n  R   contains  at  least  the  following  point 


E.  -A. 

z  =  £  .  [  ^ --1  J  e  .  +  ;,  ■     ;  e  . 

1  £  A  -£  .  ,  A  .   i  ]fi  i 

j  e  A  j   J  J    J 


n 


Usina  this  lemma,  a  criterion  for  determining  S  -stabili- 
zability  for  hyperplanes  can  be  easily  formulated. 


Theorem  2.2.2.   Given  a  hyperplane  H  =  {(a,,a~, 
An  +  ^-=iA-a-  =  0},  let  h(x)  be  the  polynomial 


,a    ) 
'    n 


n    w-n-l 


("Jx 


Let    A 


n       A     ,  n  - 1 .     a  - 1     .  nio 

.     „A . (  .     . ) aJ         i    =    0,1,2,.  .  .  ,  n. 

3=0    33-1 


j  =  l    n-j  vn-j 
Then,     the    following    statements    are    equivalent. 

(i)       H    is    RSn-stabilizable. 

ex 

(ii)   H  is  Sn-stabilizable. 

a 

(iii)   At  least  two  of  the  numbers  A.,  i  =  0,l,2,...,n,  are 
of  strictly  opposite  signs. 


(iv)       Thr-ro    exist    nomioq.it  i  ve    i  nt  ccjer.s     j    '    k    such    lli.it 

th(j)  (-a)  ](hCk)  (-a)  J  <  0  whore  h  (  ' }  (-a)  and 

(k) 

h    (-«)  are  respectively,  the  jth  and  kth  deriva- 
tive of  h(x)  at  x  =  -a. 

Proof:   (i)  °  (ii).   Recall  that  S   and  RS   are  open  connected 

sets  in  R  .   Hence,  by  lernma  2.1.2,  H  n  S   *  §    if  and  only  if 

H  n  co(Sn)  *  (1  .       Also,  H  n  RSn  x    ±    if  and  only  if  H  n  co(RSn)  x   <J>. 
a  'a  J  a 

Since  co(RS  )  =  co(S  ) ,  one  concludes  that  the  Dronerty 
a         a 

II  n  Sn  *    d)  is  eouivalent  to  H  n  RS  *    d>. 

a.        H  a 

(ii)  °  (iii).   By  lemma  2.1.2,  we  know  that  the  property 

H  n  S   *  <k>  is  eouivalent  to 
a         2 

H  n  co(Sn)  *    c|)  (2.16) 

Recall  from  (2.14)  that, 

,  „ n ,    r  r n   r  /  n  >  ,   >  j  ,  - n   ,  n -  i  >  ,   .  j  -  i   ,  -*  i    ^  „  ■, 
co(Sa)  =  {Zj  =  1[(j)  (-«)J  +  £i=i(j_i)(-a)J   r.DejIr.  >  0} 

Thus,  (2.16)  is  true  if  and  only  if  there  exist  r.  >  0, 
i  =  l,2,...,n  such  that 

An  +  Z    [(")(-a)j  +   E  (""})  (-a)J_ir  ]A   =  [A   +  Z     (n)(-a)^l 
U    j=l   ]  i=l  :  X  ]  j=l  J 

n    n  ■  .  n 

+  I    [  Z     (   X)(-a)D   A.]r.  =  An  +   Z  A.r.  =  0         (2.17) 
i»l  j  =  l  -1"1  J       L     °    i  =  l  :  X 

By  lemma  2.2.1,  the  existence  of  such  r. 's  which  satisfy 
(2.17)  is  equivalent  to  the  property  that  at  least  two  of  the 
numbers,  A.,  j  =  0  ,  .1  ,  2  ,  .  .  .  ,  n  are  of  strictly  opposite  signs. 


(iii)  -•-  (iv).   This  equivalence  follows  readily  from 
the  observation 

h(l)  (-a)  =  "fA.   for  each  i=0,l,2,...,n     [J 

ill  'III  LJ 

A  special  case  worth  mentioning  is  when  a  =  0. 

Corollary  2.2.3.   A  hyperplane  H  =  {(a  ,  a  _,..., a  )  e  Rn | 

A0  +  ^i  =  lAiai  =  °^   1S  s0~stabllizai:)le  if  and  only  if  A. A   <  0 
for  some  integers  j  *  k  where  0  <  j,  k  <  n. 

Proof:   The  proof  follows  readily  from  the  fact  that  A.  =  A. 
11 

for  all  i,  when  ex  =  0.  [] 

It  is  interesting  to  note  that  the  equivalence  between 
(i)  and  (ii)  of  theorem  2.2.2  means  that  on  any  Sn-stabi lizablc 
hyperplane  not  only  can  one  find  points  (a-,  ,  a9  ,  .  .  .  ,  a  )  e  H 
such  that  x   +  E._,a.x   J  =  0  has  only  roots  with  real  parts 
smaller  than  a  but  one  also  can  find  points  ( a  -,  ,  a-, ,  .  .  .  ,  a  )  £  H 
for  which  all  the  roots  of  the  polynomial  equation 
x   +  E._,a.x   -*  =  0  are  real,  distinct  and  less  than  a.   In 
fact,  theorem  2.2.2  can  be  generalized  to  the  following: 

Theorem  2.2.4.   Let  Y    be  a  set  of  complex  numbers  satisfying 
the  condition 


{ x |  x  <  a }  £  f  c  { z  |  Re ( z )  <  a } 


Let  n(D  =  { (a,,a2, . . . ,a  )  e    R 


0  implies 
x  e  f}.   Then,  the  property  H  n  II  (D  *  <£  is  equivalent  to 


n  ,  _  n      n  - 1 
>:   +  L  .    ,  a  .  x   J 
1  =  1  1 


H  n  RS  *    i. 
0. 


Proof:       Observe    that    I<s"    c    ii(i').       Therefore,    H    n    l->s"    *    * 

implies    H    n    11(1')     x    ,\, .       Conversely,     if    Ii    n    J|(r)     *     j, ,     then 

H    n    co(Il(D)     *    (j).       Since,    by    theorem    2.1.10,    co(Jl(r))     =    co(RSn) 

a  ' 

it  follows  that  H  n  co(Rs[J)  *  <J> .   This  implies  that 
H  n  RSa  *  d,,  because  RS    is  open  connected  in  Rn . 

2.3   Bounds  of  Stabili  zabi  1  i  ty 

Let  H  =  {(al'a2 'an)  c  Rn|  A0  +  -n=1A-a-  =  0}   be  a 

hyperplane  in  Rn .   Since  S^    c  s"  if  and  only  if  a  <  S,  it  is 

obvious  that  there  exists  always  an  a.  such  that  H  n  Sn  *  d. 

U  a 

for  all  a  >  aQ.       Hence,  the  set  {a|  H  n  Sn  ^]  is  always 
unbounded  from  above.   The  interesting  question  is  "what  is 
the  inf{a|  H  n  s"  *    <}>}?"    The  next  theorem  gives  the  answer. 

Theorem  2.3.1.   Given  a  hyperplane  H  -  {(a,, a  ,...  a  )  e  Rn I 

12'      n       ' 

AQ  +  ^=lAJaj  =  °}   in  R"'  let  hU)  =  Ej=0An-j(n-j)xn"J  and 
k  =  max{j|  A.  *  o}.   Also,  let  n  =  -max{a|  IIk   h(i)(a)  =  0). 

Then,  inf{a|  H  n  Sn  *  * }  =  n. 
a 

Proof:   Observe  that  h(x)  and  its  derivatives  h^'(x),  j  =  1, 
2,..., n, have  the  same  leading  coefficient  A     Without  loss  of 
generality,  we  assume  that  A   >  0.   Then,  for  every  a  <  n , 
h    (-a)  >  0  for  all  i  =  0,l,2,...,n.   In  view  of  theorem  2.2.2, 
this  implies  that  II  n  Sn  =  <f  for  every  n    <    i\  .       Hence, 

inf  {a  |  H  n  sJJ  *  (J.)  >  n  (2.18) 

Conversely,  by  the  definition  of  n ,  there  exists  some  positive 
integer  i  <  k-1  such  that  h(l)(-n)  =  0.   Also,  since  h(l)(x) 


is  a  polynojrun]  of  [positive  degree  with  positive  lending 
coefficient,  there  exists  an  c  >  0  small  enough  such  that  for 
every  a  e  (r,,n  +  c)  either      h(l  +  1)  (-a)  >  0   or      h  ( i  +  1)  (-a)  <0. 
In  the  latter  case,  since  h(k)(-a)  =  (nl/k!)A   >  0  (we  assumed 
Ak  >  0),  one  has  Ch (l+1) (-a) 1 [ h (k) (-a) 1  <  0  for  all  a  c     (n,n+e). 
In  the  first  case,  one  has  h(l  +  1)(-u)  >  0  for  all  a  e  (n,n+e) 
and  h  X    (-a)  =  0.   Consequently,  h(l)(-a)  <  0  for  all 
a  £  (n,n+e).   Therefore,  h(k)(-a)  h(l)(-a)  <  0  for  all 
a  £  (n,n+e).   By  theorem  2.2.2,  both  cases  imply  that 
H  n  Sn  *  $>    for  all  a  e  (n,n+e). 

Using  (2.18),  one  concludes  that  inf{ct|  K  n  Sn  *  ij> }  =  n.D 
With  n  defined  as  above,  the  assertion  of  theorem  2.3.1  can  be 

restated  as  follows:   "A  hyperplane  H  is  S  -stabilizable  if 

a 

and  only  if  a  >  n . " 

Remark:   Recall  that  in  the  introduction,  we  noticed  that  if 

m  =  n,  then  for  any  a  £  R  there  exist  always  parameters  k., 

1 

i  =  1,2,  ...  ,n,  such  that  all  the  roots  of  the  equation 
P(x)  +  z"=1kiPi<x)  =  0  (*)  satisfy  the  inequality  Re(x)  <  a. 
In  other  words,  one  may  "push"  all  the  roots  of  (*)  toward 
the  left  side  of  the  complex  plane  as  far  as  one  wishes  by 
the  proper  choice  of  k.'s.   The  theorem  above  shows  that  this 
is  not  true  for  the  case  m  =  n- 1 .   In  fact,  the  theorem 
describes  the  "limit"  which  can  be  achieved  for  all  the 
possible  choices  of  k. 's. 

Before  we  prove  the  next  theorem,  let  us  divert  for  a 
moment  to  prove  the  following  property  of  polynomials. 


I  i 

Lemma  2.3.2.   Suppose  that  A  is  the  only  one  (counting  the 
multiplicity)  positive  real  root  of  a  polynomial  f  whose 
derivative  f  has  at  most  one  positive  real  root.   Then, 
f *  (x)  *  0  for  all  x  >  A . 


Proof:   Assume  the  contrary:   there  exists  xQ  >  A  such  that 
f(xQ)  =  0.   Then,  either  f " (xQ)  =  0  or  f " (xQ)  >  0  or  f " (xQ) <  0 
We  will  show  in  the  following  that  all  these  cases  load  to 
contradictions.   Without  loss  of  generality,  we  assume  that  f 
has  a  positive  loading  coefficient. 

(i)   If  f"(x0)  =  0,  then,  by  the  assumption  f  (xQ)=0rone 
finds  that  x  -  xQ  is  a  positive  double  root  of  f'(x)  =  0. 
This  contradicts  the  hypothesis  that  f  has  at  most  one  posi- 
tive real  root. 

Note  that  A  is  the  only  one  positive  real  root  of  f(x) 
and  that  f(x)  has  positive  leading  coefficient.   Hence, 

for  all  x  >  A.   Therefore,  there  exists  c  >  0  small 

and  for  all  x  e  ( A, A  +  e) 


f(x)  >  0 

>gi 


enough  such  that  A  +  €    <  xQ 


f ' (x)  >  0 


(ii)   Suppose  that  f " (xQ)  >  0.   Then,  from  the  assumption 
f" (x  )  =  0,  it  follows  that  there  existsa  6  >  0  small  enough 


such  that  A  +  r. 


and  f ' (x)  <  0  for  all  x 


x0-6,xQ) 


(A  +  -)  >  0,  and  f'(x0-6/2)  <  0,  by  the  intermediate 


Since  f '  \  ■  =p  "    u  '  cl""  *~     v"0 


where  A  <  A  +  i  /2 


value  theorem,  there  exists  x-^  where 

such  that  f(x1)  =  0.   Therefore,  f  contains  two  positive 

real  roots  x  =  x1    and  x  =  xQ,  a  contradiction. 


6/2 


(iii)   Suppose  that  f" (xQ)  <  0.   Then,  since  f'(xQ)  =  0, 
there  exists  a  f.  >  0  such  that  f'(x)  <  0  for  all  x  c  (x  ,x  +6 
Observe  that  f  has  positive  leading  coefficient.   Hence, 
there  exists  6  >  xQ  +  6  such  that  f'(x)  >  0  for  all  x  >  S. 
Consequently,  f (xQ  +  6/2) -f  (6)  <  0.   By  the  intermediate 
value  theorem,  there  exists   x0  where  xn  <  >•   +  6/2  <  x   <  R 
such  that  f'(x2)  =  0,  a  contradiction.   This  completes  the 
proof.  n 


Theorem  2.3.3.   Given  a  h 


yperplane  H  =  {  (a. , a_  ,  .  .  .  ,a  )  e  Rn 
1   /      n 


A0  +  E.  jA.a.  =0},  let  k  =  max{j|  A,  *  0}.   Let  h(x)  = 
^"=0Aj(")xD.   Let  n  =  inf{a|  H  n  s"  *  <J>}.   If  there  exists  a 
nonnegative  integer  P..  <  k-1  such  that  A.  <  0  for  all  j  <  £ 
and  A.  >  0  for  all  j  >  I ,    then  n  =  -max{a|  h(a)  =  0}. 

Proof:   Obviously,  -max{a|  h(a)  =  0}  >  -max{a|  nk_  h ^ *  (a)  = 0}=n. 
Conversely,  recall  that  h(x)  =  7n      A.  [U)y?  .       The  hypothesis 
implies  that  there  is  exactly  one  variation  of  signs  among 
the  coefficients  of  h  ^  (x)  for  each  0  <  j  <  £-1  and  no 
variation  of  signs  of  coefficients  of  h  -  (x)  for  j  >  I. 
Therefore,  by  Descartes  Rule,  the  polynomial  h  -1   has  exactly 
one  positive  real  root  if  j  <  2.-1  and  no  positive  real  root 
if  j  >  I.       It  follows  that 

£  - 1  k     . 

-max{a|  li   h(^(a)  -  0}  -  -max{a|  i;   h(;i)(0  =  0}  =  n 
1=0  j^0 

Furthermore,  by  applying  lemma  2.3.2,  to  h,h';  h',h ",..., 

(9-2)      (  £-1) 
h     ,k  '    ,  one  obtains  that,  for  each  0  <  j  <  £-1, 


maxlcxl    h{j)  (a)     =    0}    <    max{a|    h(^    i: 


0} 


Consequently , 
f-1 


(J) 


max{a|  II   hvj;(a)  =  0}  =  max{a|  h(a)  =  0} 
j=0 


Hence , 


H  =  -max {a |  h  (a)  =  0} 


Remark:   Given  a  hyperplane  H  c  Rn ,  suppose  -inf{a|  H  n  Sn  x   4 ,)  =n 
From  the  point  of  view   of  applications  it  is  important  to 
actually  find  the  points  in  H  n  sj|  <  <f>  if  a  >  r,  .   Unfortunately, 
we  do  not  have  a  general  method  for  such  a  task.   Following  is  a 
special  case  for  which  we  do  have  a  method.   Suppose  H  is  a 
hyperplane  in  R1  such  that  n  =  -max{a|  h(a)  =0}.   Then, 
h(-n)  =  0.   Recall  that  h(x)  =  En=QA.(n)xj.   Therefore, 
Ej-qAj  (!?)  (-n)3  -  0.   This  implies  that 

((")  (-n),  (n2)  (-n)2,...,(^)  (-n)n)  £  H 

Note  also  that 


((?)  <-n),  C^)  <-n)  2 (n)  (-n)n)  e  rs" 

1  ^  n  n 


n    n  n    n 
identical  roots  x  =  n.   Hence, 


since  the  polynomial  x"  '  E"  , (") (-n) 3xn  J    =  0  has  only 


n 


((  )  (-M),  (")  (-n)" (")  (-r,)n)  e    H  n  RSn 

i-                      /■                               n  n 

We  remark  that  this  method  can  be  applied  to  the  kind  of 
hyperplane  described  in  theorem  2.3.3. 


In  the  following,  we  illustrate  the  results  which  we 
obtained  in  this  chapter  by  an  example  chosen  from  Anderson 
et   al.  [1]. 


Example:   Consider  a  system 


E)   x  = 


1 


0  1  0 
0  0  1 
0   13    0 


f   1 
0 

0 

1 


y  = 


0    5-1 
-1  -1    0 


x  ,    u  =  Ky  =  [  v  ,  w  ]  y 
Then,  the  characteristic  polynomial  of  the  closed  loop  system 


s3  +  vs2  +  (w-5v-13)s  +  w  =  0 


2.19: 


Let  II  denote  the  set  of  all  the  coefficients   vectors  arise 
from  all  possible  choice  of  K  =  [v,w],  i.e., 


H  =  {(a,,a„,aJ|  a,  =  v,  a?  =  w-5v-13,  a-.  =  w,  for  some 


v,w  e  R }   =  {(a,  ,a~,a.j|  5a,  +  a . 


+  13  =  0 } 


Note  that  H  forms  a  hyperplane  in  R  .   Since  A,  =  5  >  0  and 

A   =  -1  <  0,  by  corollary  2.2.3,  the  system  (E)  is  SQ-stabili- 

1x2 

zable.   In  other  words,  there  exists  a  K  e  R     such  that  all 

the  roots  of  (2.19)  have  negative  real  parts. 

Consider  the  polynomial  h(x)  defined  in  theorem  2.3.3. 

3      2 
For  system  (E),  h(x)  =  -x   +  3x   +  15x  +  13.   By  theorem  2.3.3, 

one  concludes  that 


in f  {a  I  H  n  S3}  =  -max ht  |  h(ct)  =  0 


■5.91 


This  moans  that  no  matter  how  wo  choose  K    [v,wl,  t.horo  exists 
at  least  one  root  of  (2.19)  with  real  part   qreater  than  or 
equal  to  the  number  -max {a  c  R|  h(a)  =  0}  ss  -5.91. 


'  a ,  8 


CHAPTER  III 
-STAB  I  LIZ  ABLE  11 YPERP  LANES 


Given  a  <  6,  let  C   ,  denote  the  set  of  all  the  points 

Ot  ,  D  l 

(a, f a2 , - . . , a  )  e  R   for  which  the  polynomial  equation 

x   +  ^i  =  1  aix       °  lias  only  roots  inside  the  circle  centered 

at  ((a+8)/2,0)  having  radius  (6-a)/2.   Let  RCn  _  be  the  set 

a ,  p 


i-  ( a -,  /  a  ,  .  .  .  ,  a  , 
12       n 


e  Rn|  all  the  roots  of  xn  +  Ln    ,  a . xn  J    =    0 

3  =  1   J 


are  real,  distinct  and  strictly  betwen  a  and  8}. 


A  hyperplane  H  c  R   is  said  to  be  stabilizable  with  respect 

to  the  circle  {z|  |z  -  (B+a)/2|  =  (6-a)/2)  if  H  n  Cn  „  *  <t> . 

a ,  b 

H  is  stabilizable  with  respect  to  the  interval  {x|a<x<3}  if 
H  n  RC   -  *  (J> .   In  this  chapter  we  investigate  the  following 
two  problems: 

(i)   What  are  the  co(RCn  0)  and  co(Cn  „)? 
(ii)   How  can  one  characterize  the  hyperplanes  which  are 
stabilizable  with  respect  to  a  given  circle  and  the  hyperplanes 
which  are  stabilizable  with  respect  to  a  given  interval? 

Particularly  interesting  from  the  stability  point  of  view, 
is  the  cise  of  the  unit  circle  centered  at  the  origin.   In 
fact,  the  question  of  the  stabiliznbili ty  of  a  hyperplane 
with  respect  to  the  unit  circle  centered  at  the  origin  arises 
from  the  problem  of  direct  output  feedback  stabilization  of 
a  sampled  data  control  system.   The  well  known  Schur-Cohen 
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criterion  which  characterizes  C_,  .  consists  of  a  system  of 

nonlinear  inequalities.   Hence,  it  may  be  a  difficult  task 

to  answer  questions  (i)  and  (ii)  directly  in  terms  of  those 

inequalities.   Instead,  by  an  invertible  linear  transformation 

constructed  in  section  1,  we  shall  reduce  (i)  and  (ii)  to 

problems  solved  in  chapter  2.   As  a  main  result,  we  prove 

that  co(C   „)  and  co ( RC   „ )  are  both  equal  to  the  interior  of 
a,  8  a  ,  B 

the  simplex  generated  by  the  points  (q .,, q .„,...,  q  .  )  e  R 

i  =  0,l,2,...,n  where  q..  are  the  coefficients  of  the  polyno- 


mial  (x-a)    (x-8)   =  x   +  I  q^x 


n-i 


3.1   The  Convex  Hulls,  co(C  Q)     and  co(RC  0) 
We  will  denote  the  standard  basis  of  R   by  { e . } .  ,  . 
For  brevity,  if  there  is  no  confusion  of  the  dimension,  we 
will  simply  write  {e.}.  ,  .   First,  let  us  extend  lemma  2.1.5 
to  the  following  sets 


n     ->    +1  n 

S  =  {  Z   a.e"    I  a   >  0,    I       a.x 

j=0   3  3+1'   0        j=0   ^ 


n-i 


0  implies  Re(x)  <  0} 


n       +  1  n 

RS  =  {  I      a.e.  .  I  a„  >  0,  all  the  roots  x  of    7. 


n-i 
a.x   J  =  o 


are  real,  distinct  and  x  <  0} 

Lemma  3.1.1 .   co  ( S )  =  co ( RS )  -  r"  +  * . 

Proof:   Obviously,  by  Ilurwitz  criterion  and  the  convexity 
property  of  r"  +    one  has  co(RS)  c  co(S)  c  R+ 


.  _  „n      ->-n  +  l    nn+l   .  , 
Conversely,  if  E .  n    a .e . , ,  e    R    ,  then 
J        -|=0   ]  3  +  1     + 


j  =  l 


(a  ./a.) e . 
D   0  j 


Recall  that  r"  =  co(Rs").   Therefore,  En   (aVan)e"   can  be 
+         0  3=1  J       u   j 


expressed    as    a    convex    coinbinat  ion    of    [joints    from    RS    .       In 


S=i(biJ/ao)ej 


RSQ  and 


other  words,  there  exist  point 

parameters  k-  >  0  satisfying  E  .  k  .  ■-  1,  such  that 

E  .  ,  a./an  e.  =  Z  .     , L E .  .  k.(b../art)je.  or  equivalently, 
j  =  l   j   0   j     j  =  li  =  liijO    j  J 

■n  +  1 


rn      +n+l      ->-n  +  l  ,  rn   rvn   ,  , 

E  .  A  a  .  e  .  ,  ,  =  a  n  e .    +  E  .  ,  [  E  .  .  k  .  b  .  .  e  .  ,  , 
3=0   3  3  +  1     0  1       j=l   1=1  l  i]   j  +  1 


Observe  that  for  each  i,  the  equations  x   +  E.  ,  (b.  -/a„)x 

i=i  n'  o 


(3.1) 


bj^x     —    0  have  the  same  roots.   Hence,  from 


and  a„x   +E.  ,  ,  ■ .  . 
0       :=1   13 

rn   /u   /   \  ^n    n^n        1  i.        -»*n+l  ,  vn    .        »-n+l     ~ 
E  .  ,  b.  Van  e.  e  RS- .  one  obtains  ane,    +  L .  n    b .  . e .  , ,  e  RS 
]=1   13'  0   ]      0'  0  1       j  =  0   13  3+1 

Therefore,  by  (3.1),  one  concludes  that  E.  , a  •  e  .  , ,  e  co(RS) 

3=1  3  3+1 

for  any  E.  ^ae.,-,     e  R,   .   This  implies  R,    c  co(RS).   Q 
3=0  3  3+1     +  +    -    v   ' 

For  any  fixed  a  <  8,  define 

n  +  1 


C   _  =  IE .  n  a .e . , 

a ,  0      ]  =  0   3  3+1 


r.n  +  1 1 
R      a 


, .  n       n  -  1 
E  .  n  a  .  x   J 
1  =  1   3 


0  implies 


(a+B)/2|  <  (3-a)/2) 


and 

RC   „  =  IE.  „  a.e.,,  e  R     a,,  >  0,  all  the  roots  of 

a, 6      3  =  °   3  3  +  1  ° 


3=0  "j 

a  and 


=0  are  real,  distinct,  and  striclty  between 


Let  T   „  be  the  transformation  on  R 

a ,  6 

rn     ->-n  +  l    nn  +  l 
E  .  n  a  •  e  .  .  ,  r  R 
3=0   3  3+3 


n  +  1 


such  that  for  each 


T   _  (  ),  a.e.,,)  =    L  a  .  c .  . 


where    a.,     i    =    0,l,2,...,n    are    obtained    from    the    following 
3 

polynomial    expansion: 


J!j 


I    a.  (s-a) 
j  =  0  J 


n""D(s-B)]  =  I    a.sn"j 

J-0  3 


For  a  =  1  and  8=1,  it  is  well  known  that  T  .    is  a 

""■*■*  J- 

linear  transformation  which  maps  S  onto  c",    (Hermite  [13, 
Marden  [1]).   The  matrix  representation  Q      of  T     ,  due 
to  its  application  to  the  stability  of  a  discrete  system, 
has  been  investigated  by  many  researchers  such  as  Power  [1], 
Fielder  [1],  Halijak  and  Moe  [1],  Jury  and  Chan  [1],  Duff in  [1], 
etc.   Barnett  [1]  has  studied  the  linear  transformation  T 
induced  by  a  more  general  bilinear  transformation  z  =  j^y 
and  showed  that  the  entries  of  the  matrix  representation  of 
T  can  be  obtained  by  an  interesting  algorithm.   For  our  purpose, 
we  need  the  following  special  case  of  Harnett's  (i.e.,  6=1 
and  y  =  1) • 


Lemma  3.1.2.   For  any  a  <  3,  T     is  an  invertible  linear 

transformation  on  R     which  satisfies  T   , (S)  =  C   ,  an 

a  ,  6        a  ,  6 

T   .  (RS)  =  RC"  0  . 

Proof :   We  divide  the  proof  into  three  parts: 

(i)   T   n  is  linear, 
a ,  d 

„       rn     ->n  +  l    ,  yn      ,  -+n+l   ,    rn   ~  ►n+1 
For  any  I.    na.e.,,  and  l.       b.e.,,,  let  l.     .c.e,, 
]=0  j    ]+l       n=0  ]  ]+l'       n=0  ]  3+1 

T  c(ln       (a  .fb.)entJ)  .   By  the  definition  of  T   ,  ,  the  c., 
C/B   1=0   J  j       n+1      J  a, 6        ] 

j  =  0,1,..., n,  are  obtained  from  the  following  polynomial 
expansion 


j=0 


•:  c  -sn  D  =    l     (a.+b.)  (s-a)n  :(s-6):  (3.2) 


j  =  0 


T,   rn    ~  -►n+l     ,   vn    ~  vn  +  l 
If   E  .    n    a  . e  .  .  ,   and   E  .     n    b 


-n+1 


_,  .  e  .  ,  ,  arc  the  images  of  E.  „  a.e.,., 
]=0   ]  3+1         3  =  0   33+I  J        3=0   ]  3+1 

and  I  ,    „  b.e.,,  under  T   „,  respectively,  then  the  right  hand 
j  =  0   j  j  + 1         a ,  B 

side  of  (3.2)  equals  En  =  Q  (a  .  +  b, )  s11-^  .   Thus,  c   =  a.  +  b  .  , 

j  =  0,1,2, . . . , n .   Hence, 


n  n         ,  n        , 

T  A    Z    (a.+b.)en;J)  =  T  A    I      a.  en;!)  +  T  A    I      b.e"*,) 

a,Bvj  =  Q   3   3   1  +  1      <*,B  ..=n   3   D  +  -!      a>B  ^-n   '  n  '  > 


j  =  0 


3  =  0 


(ii)   T     is  invertible. 
ex ,  p 

To  show  that  T   .  is  invertible,  since  T     is  linear, 
a,  6  a, 6 

it  suffices  to  show  that  T  n    is  onto.   In  other  words,  for 

a,  B 


-n  +  1 


any  E.  _  a.e.,,  we  need  to  find  E.  n  a.e.,, 
1   3=0   3  3+1  3=0   3  J+1 


n+l        U       *-l    u 

a.e.,-,  such  that 


I   a  .  (s-a)n  j  (s-6) j  =   E 
j=0   D  j=0 


n-3 


(3.3; 


Substitute  (6z-a)/(z-l)  for  s  in  (3.3)  and  multiply  both  sides 
by  (z-l)n.   Then,  (3.3)  becomes 


E   a . zn  ^    =   E   a  .  (Bz-a 


n-3 


z-1) 


3.4; 


3  =  0 


j  =  0 


In  view  of  this,  one  may  simply  choose  a.'s  to  be  the  coeffi- 
cients of  the  polynomial  expansion  of  the  right  hand  side  of 

(3.4).   This  proves  that  T  .    is  onto.   Hence,  T     is 

a ,  B  a  ,  p 


invertible. 

(iii)       T  (S)    =    C 

a  ,  B  ' 


xfB      and   Ta,P(RS)    =    RCa,6 


Consider    the    1-1    transformation 


s-P 


(3.5) 


which  maps  the  interior  of  the  circle  {s|  |s-(n+fl)/2|  =  (B-u)/2) 
onto  the  open  left  half  plane  {z  |  Rc(z)  <  0}.   In  particular, 

(3.5)  maps  the  interval  {s|a  <  s   <  6)  onto  the  negative  real 
axis.   Hence,  for  any  nth  degree  polynomial  p(z) ,  the  state- 
ment: "all  the  roots  of  p(z)  =  0  are  in  {z  |  Re(z)  <  0}"   is 
equivalent  to  the  statement:   "all  the  root?  of  p(s)  = 

(s-6)np(S4|)  =  0  are  in  {s|  |  s-  (a+B)  /2  |  <  (6-a)/2}."   Further- 
more, {z|  p(z)  =  0}  c  (z|Re(z)  <  0}   if  and  only  if 

{s  |  p(s)  =  0}  c  {s|  a  <  s  <  6).   Observe  that  for  any  p(z)  = 

n  .      .     n 


Z  .  „  a  .  x 
D  =  0   D 


n-i       ,       D  s  n  ,s-ti,  _ 


j=0 


a  .  (  s  -  a ) 


n-j 


s-B 


l      a  .  s 


n-j 


j  =  0 


Recall  that  this  was  how  the  image  of  £•  «  a^en  +  i   under  T^  ~ 

was  defined.   Hence,  T  ol    determines  a  one  to  one  correspon- 

a ,  p 

dence  between  the  Hurwitz  polynomials  and  polynomials  whose 

roots  lie  inside  the  circle  {s|  |s-(a+p)/2|   =  (£-cx)/2). 

Similarly,  T   ,  determines  a  1.-1  correspondence  between  the 
1        a ,  p 

polynomials  which  have  only  distinct  negative  real  roots  and 
the  polynomials  of  which  all  the  roots  are  distinct  real  and 
strictly  between  a  and  B-   If  the  leading  coefficient  aQ  of  a 
Hurwitz  polynomial  is  negative,  then  a.  <  0,  j  =  1,2, ...,n. 
Hence,  aQ  =  z"  =  0  a.  <  0  which  means  that  Ta/ g (^ j=oaje j+1^  = 


^0^e^}  *  C^f.   Therefore,  ^(S)  =  co(C^  R)  and 


fX  r 


T     (RS)  =  Re" 


a,P 


a 


Following  a  similar  technique  as  Barnott  I  11,  we  compute 

ex,  6 


explicitly  the  matrix  representation  Q^    of  Tq  _  relative 

.n+1 


to  the  usual  standard  basis  of  R 

Note  that  Q   ,  is  an  n+1  by  n+1  matrix  such  that 

ex ,  D 

Vb'JoVj  +  I1    =    <Val anIOc.,f 


ih 


Hence,  from  the  definition  of  T   .,  for  cacti  fixed  i,  the  row 

CX  ,  D 

^qiO'qil' ' - * ,qin)    of    Qa    8      roPrescnts    the    coefficients    of    the 
polynomial 


,         ,  n-i  ,       .,i 
(s-a)  (s-6) 


(3.6! 


By  a  simple  combinatorial  argument,  one  may  expand  (3.6)  as 

(s-a)n-i{s-6)i  =  I    C(-l)j   ^  (i)  (n~j)a^-k6k]sn-^    (3.7) 
j=0        k  =  0  k   D  K 

Hence,  for  each  0  <  i  <  n  and  0  <  j  <  n,  one  obtains 

q-,  =  (-D3       £  (b   ("~£)aj'V  (3.8) 

ID  k=0  K   D~K 

In  particular,  from  (3.7),  it  is  easy  to  see  that  q,n    =  1, 

lO 

i  =  0,1,2, ... ,n. 

We  remark  that  the  determinant  of  Q  a    is  equal  to 
n(n-H)  a'b 

2 
(a-B)        (see  appendix  II)  . 

Now  consider  a  transformation  T   „  on  R     defined  by 

a  ,  8  x 

n      ,i      n      ,  , 

T  AT,    S.en!h  =   E  a.en!} 
a,B  j=Q  ]  3+1     j=Q  1  1+1 

where  the  a.,  j  =  0,1,2, ...,n  are  obtained  from  the  following 

polynomial  expansion: 

n  _  .      .     n 

T.    a.(Bz-a)n  D(z-1)D  =  Y.    a.zn_:> 
j=0  ]  j=0  ] 

6  7  ~  o. 
By  the  transformation  z  — >  s  =   J_-,   and  with  the  same 

technique  used  in  establishing   lemma  3.1.2,  one  can  show 

that  T  0    is  a  linear  invertible  transformation  on  R     which 
o. ,  P 


jy 


satisfies  the  relations  T  (  < ;        )    =    rs  and  T   . (cn  „)  =  S 

Furthermore,  one  finds  that  the  entries  q. .  of  the  matrix 

representation  0„  D  of  T  „  are 

CX  ,  p         CX  ,  p 


ID 


k=o  k  j~k' 


We  remark  that  one  may  use  the  algorithm  proposed  by  Barnett  [1] 

to  generate  the  entries  of  matrix  Q 

a  ,  (3 

For  the  unit  circle  case,  i.e.,  a  =  -1  and  6  =  1,  it  is 
interesting  to  note  that  the  matrices  Q  ,    and  Q      which 
are  associated  with  T_,  ^  and  T  ,  -,  ,  respectively,  are  equal. 
Duff in  [1J  and  Jury  and  Chan  [1]  have  computed  the  matrix 
explicitly  and  have  proved  the  following  relations  between  the 
entries  of  Q_-,  ■, 


qi3  =  qi-i,j  "  qi-i,j-i  "  qi,j-i  x  ■  1'   i   s  L 

Q-1,1Q-1,1  =  ("2)n(n+1)I 


o.io: 


Obviously,  (3.10)  can  be  easily  programmed  on  a  computer  to 
generate  the  matrix  Q_  -,  -,  . 

Now,  recall  that  if  T  is  an  invertible  affine  transfor- 
mation on  R  ,  then  for  any  set  S  c  R  ,  co(T(S))  =  co(T(S)). 
Using  this,  we  obtain  the  following  proposition. 

Proposition  3.1.3.   For  any  fixed  a  <  8  and  every  n,  the  sets 
co(C\,  „)  and  co  ( RC   „)  arc  both  equal  to  the  polygonal  convex 


a,  h 


,n  +  l 


set  in  R  formed  by  the  intersection  of  the  following  half 

spaces : 

nT-[(-DJ  E  ("-j)(j)a^en-^+k|a.  >0   j  =  0,1,2 n. 

i=0  k=0  '  1 


4  0 


Proof:   Recall  that  T   „ (C 
a ,  p   a , 


=  S  and  T   .(RC"  Q. 


From  lemma  3.1.1,  it  follows  that 


RS 


T   „(C0(C"  „))  =  co(T„  H(C"  B))  =  CO(S)  =  Rn+1 


a, B     a, 6 
also, 


ct,B   n,6' 


1   „(co(Rc"  p))  =  co(Ta^(RC^e))  =  co(RS)  =  R"  +  1 


a,B 


'a, 6' 


-n  +  1 


Hence,  L^=Q    ajej+1  e  co(Ca^R)  =  co(RCa  P)  if  and  only  if 
Ta,S(Lj=0  aje]+l)  =  (a0'al an)Q«,6  >    °  °r  equivalently , 


i=0L  (  l)     ^k=Olj-k; (k,a    P 
i  =  0,1,2, ...,n. 


ja.  >  0  for  each 


Observe  that,  in  terms  of  the  usual  basis  vector  of 

Rn    the  positive  orthant  R^+   can  be  expressed  as 
-n  +  1 


{Zi=0  ki^i+l  I  ki  >  0)-   Xt  follows  that 


co(c"  p)  =  co(Trv  p(S))  =  T   a(co(S))  =  T   D(R"+1) 


'a,B 


ri, 


a ,  B  v  + 


=  {  E  k,  T   D(e?!J")  !  k,  >  0} 


i  =  0 


'i   a,Bv  i  +  1 


{  Z  k.(q    q.1# 
i  =  0 


-,qin)|  k.  >  0} 


where  (q .Q ,q.. j , . . . , q.  )  is  the  ith  row  (0  <    i  •  n)  of  matrix 
Q   0.   This  gives  us  another  characterization  of  co(Cn  ., )  and 

co(RCn  .) . 

a,  B 


Proposition  3.1.4.   For  any  fixed  a  <  B  and  n,  co(Cn  „)  and 

o. ,  p 

co(RC  Q)    are  equal  to  the  interior  of  the  convex  cono 
oc ,  fc> 


41 

■ir       ,  ,  ,  •  n  Mil  1         ,  n  +  1       .  .  ,  • 

generated    by    the    vector::    i.  .    n    q.  ,c.     ,     '     R  with    oriuin    as 

-' "  :]-0        i]     j  ll 

vertex,  where  for  each  i,  the  q—,  j  =  0,1,...,  n  are  the 
coefficients  of  the  polynomial  expansion  (x-cx)    (x-P)   = 

For  M  c  Rn    and  a  e  R,  let  M   denote  the  subset  of  M 

a 

defined  by  {  (aQ ,  a1 ,  .  .  .  ,  an)  c    m|  aQ  =  a}.   Note  that  Ca  [C, 

Cn  „.  RCn  -  and  RCn  „  are  related  by  the  following  relations: 
a,  B     a, 3       a, p 

(Cn  .).  =  (1)  x  c"  , 

a  ,  6  1  ex,  8 

(RCn  R) .  =  {1}  x  RCn 
a  ,  B  1  a ,  p 

The  next  lemma  shows  the  relation  between  their  convex  hulls. 
This  will  enable  us  to  answer  the  first  question  which  we 
raised  in  the  beginning  of  this  chapter. 

Lemma  3.1.5.   (i)   (co(c"  R) ) .  =  {1}  *  co(c"   )  and  (co(RC"  )), 
a,pi  a ,  p  ex  ,o   r 

{1}  x  co(RC^e)  . 

(ii)   CO(C^B)  =  co(RC^p). 

Proof:   (i)   Every  point  from  (co(C   g) ) , ,  is  of  the  form 

£k. (b. ft,b. , , . . . ,b.  ) 

.  i   lO '  ll     '  in 
l 

where  (b .  .  ,  b .,,...,  b .  )  e  c"  ,,  and  the  k.'s  satisfy  the 
1O1I       in      «,[.'  i 

relations  k.  >  0,  E.k.  =  1  and  £.k.b.n  =  1.   Since,  for  each 
1         11  i  i  i  U 

(bi0'bil bm»  (  £a,B'  bi0  "    °'  lL  follows  that 

Eki(bi0'bil bin}  =  {1}  X  ^kibi0(bil/bi0'  bi2/bi0 

bin/bi0)  c     (1)  x  co(C^B). 

The  last  inclusion  follows  from  the  definition  of  c°(C   „). 


A  2 


Conversely,  if  (  1 ,  a  ,  a  ?  ,  .  .  .  ,  a  '. 


co (C   „ ) ,  then 
ex ,  p 


(l,a, ,a_, . . . ,ar  )  =  {1}  *  { Ik . (b . , , b . 0 , . . . ,b .  )}  = 
1   z       n  l   1 1   l  Z  in 

l 

Zki(1'bil bin}  whore  (bil'bi2 bin)  c  C«,6  and 

E.k.  =  1  with  k.  >  0.   Since  for  each  i,  (l,b.,,b.  -,,...  ,b.  )  6 
ii  i  ii   i2     '  in 

C^  g;  this  implies  that  (l,a .,...,  a  )  c  (co(C   o))-,«   Therefore, 

(co(C   „)),  =  (1)  x  co(C   _) .   The  same  proof  holds  for 
a ,  8   1  a ,  B 

(co(RCn  .)).  =  {1}  x  co(RCn  .)  . 

a ,  p   1  a  ,  p 

(ii)   Since  co(Cn   )  =  co(RCn  „),  by  (i),  it  follows  that 

{1}  X  Co(Ca,B)  =  (co(£M»»l  =  (co(<,e»l  =  {1}  *  c°<<fB)- 

This  implies  that  co(Cn  D)  =  co(RCn  „).  D 

a  ,  8  a  ,  B 

Combining  the  previous  lemma  and  proposition  3.1.3,  we 

obtain  the  description  of  the  convex  hulls  of  C   D   and  RC  0. 

r  a, 8  a, 8 

Proposition    3.1.6.       Co(Cn       )     and    co(RCn    R)     are    both    equal    to 

Ua,fa? a)    e    Rn  |     (-1)  j  (")  a j  Bn" j    +      I    !(-l)j      I    (!?"*) 

n  D  i=l  k  =  0    : 

(j)aj-k6n-i-j+k3a.     >    0,    j    =    0,1,2 n} . 


Proof :   Apply  proposition  3 . 1 . 3  with  a   =  1 .  D 

Consider  the  matrix  0  n    -    (q-  •)  associated  with  T  „ 

-a,  8      ID  a'  P 

Recall  that  q.n  equals  one  for  each  i ,  0  <  i  <  n .   From 
1  i  0   J 

proposition  3.1.4,  it  follows  that 


{1}    x    co(C^B)     =    (co(C^B))1 


=    f.znki(1'qil'qi2,""qin,|ki    >    °}1 

i  =  0 

=    {1}    *    {."    ki(t3il^i2 qin,lki    >    °'     °Ci 

1  =  0  i 


1) 


Hence, 


co(C   J  =  {£.  _k  .  (q.  .  ,  q  •-,,...  ,q  •  )|k.  >  0,  Ek  . 


1} 


Proposition  3.1.7.   co(C   „)  and  co(RC   R)  are  both  equal  to 
the  open  simplex  generated  by  the  n  +  1  points  (q  ■■,,  q  ■-,,...,  q .  )e 
R  ,  i  =  0,1,2,.. . ,n  where  q. .  are  the  eoefficients  of  the 


(x-6)1  =  xn  +  Ln    .  q.  .xn  D;  explicitly, 

3=1  Jij 

D 

In  particular,  if  a  =  -1  and  6  =  1,  then  one  obtains 
the  followina  conclusions: 


polynomials  (x-a) 

/  i  \  J  v  J   * n-i wiv  j-kDk 
qiD  =  (-1)JZk=0(3-k)(k)a    6 


Corollary  3.1.8.   co(C_,   )  and  co ( RC_ ,  ,  )  are  both  equal  to 

the  open  simplex  oenerated  by  the  n  +  1  ooints  (q .  ,  ,q .  _. ,  .  .  .  ,q .  ) 

e  R  ,  i  =  0,1,2,. . . , n  where  q . .  are  the  coefficients  of  the 

i_l 

t     •  i   ,  .T.n-i,   ,  .  i     n  .  _ n       n  -  j     ,  .  .  ,  , 
polynomials  (x+1)     (x-1)   =  x   +  l  .  __ .  q  .  .  x   J  explicitly, 

"i 

q.  .  =   I     (-l)k(n~l)   (£)  . 

•^        k=0      ^k  k 

2 
Example :   Suppose  n  =  2.   Then,  co(C"   ,)  is  the  shaded 

region  in  the  figure  3.1,  and  is  bounded  by  the  following 

three  lines : 


l  +  a;L  +  a2  =  0,   1 


0 ,   1  -  a,  +  a? 


0 


In  this  case,  the  convex  hull  of  C      is  actually  equal  to 

2 
C_,  ,  itself.   Since 


"-1,1 


1     2     1 

1     0    - 1 
1-2     1 


Figure    3 . 1 


by  corollary  3.1.8,  this  implies  that  co(C_    )  is  the 

2 
triangle  in  R   having  the  points  (2,1),  (0,-1)    id  (-2,1) 

as  vertices. 

In  fact,  one  may  generalize  proposition  3.1.7  to  the 

following  theorem. 


Theorem  3.1.9.   Given  a  <  8,  let  1   D  be  the  open  interval 
■ a. ,  p 

(a, 8).   Let  V       „  be  an  open  disk  on  the  complex  plane  with  a 
and  8  as  the  end  points  of  a  diameter,  i . e .  ,  {z|  |z-(8+a)/2|  < 
(8-a)/2}.   Let  T  be  a  set  of  complex  numbers  such  that 

Ve  C~    r  £  Va,B'       Let   (F)  =  {(al'a2 an>  £  R"  I  x"  + 

E"=0a.xn  D  =  0  implies  x  e  r)  .   Then,  co(n(|'))  =  co(C,n  „)  = 

co(RCn  0)  . 
a,  8 

Proof :   The  hypothesis  on  r  implies  that 

RCa,8  C~    «{V)    £  Ca,B 
Consequently, 

co(Rc"   )  c  coU(r))  £  co(c"  R)) 
a ,  p  a ,  p 

By  the  relation  co ( RC   D)  =  co(C   D),  one  concludes  that 

a ,  8         a  ,  8 

co(tt(D)  =  co(c£  6)  =  co(c"  g).  D 


3.2   C   » - s tab i 1 i z ab 1 e  Hyperplanes 
Now,  we  consider  the  problem  of  characterizing  hyperplanes 
which  are  stabilizable  with  respect  to  a  given  circle  in 
the  complex  plane. 


4(> 


In  the  following,  for  any  hyperplano  II 


a,  ,a. ,  .  .  .  ,a. 


'-     R  I  AQ  +  I    Aiai  =  0},  we  denote  by  H  the  set  {  (a   a,  ,  .  .  .  , 


i  =  l 


.n+1 


an)  c    R    |  Z   A. a.  =  0}.   Note  that  H  is  a  hyperplane  in  Rn+1 
i=l 

Lemma  3.2.1.   Suppose  H  is  a  hyperplane  in  Rn .   Then,  for 

any  fixed  a  <  6,  the  property  H  n  c"  0    *  <j)  is  equivalent  to 

ft ,  t>  -1 

H  n  P  g  *  ♦. 


Proof:   (=»)  .   If  (a1,a2,  .  .  .  ,an)  e  H  n  CQ  _,  then  by  the 
definitions  of  H  and  C^  g  respectively ,  (1,3-,,..., a  )  c  H  and 
(l,a1,  .  .  .  ,an)  f  CQ  g.   Hence,  (1 ,  a.^  ,  a2  ,  .  .  .  ,  a  )  c  H  n  cn  „. 

^  •       If  (a0,a1,...,an)  e    H  n  C^  B  ,  then  (a^  ,a2/aQ  

an/aQ)  e    H  n  (£  g  .  D 


Using  this  lemma  and  proposition  3.1.4,  we  obtain  the 
following  criterion  for  determining  the  s tabi li zability  of  a 
hyperplane  with  respect  to  a  given  circle. 


Theorem  3.2.2.   Given  a  <  B,  let  Q     be  the  matrix  associated 
~~ a  ,  b 

with  the  transformation  T   p  (see  (3.8)).   Suppose  H  =  {(a,, a-,. 

ft  i  P  L        Z 

...,an)  e  R  |  AQ  +  £.  jA-a.  =0}   is  a  hyperplane  in  Rn .   Let 


Aif  i  =  0,1, 


,n,  be  the  numbers  obtaind  from  ( An  ,  A,  ,  .  .  .  ,  A  )   = 

0   1     '  n 


3yJ 


Qa,B(VAl V'  explicitly,  A.  =  E  j=()[  (-1)  J  Z^Q  ( »_£)  (J) 

B  ]A.  .   Then,  the  following  statements  are  equivalent. 


j-k 


(1)  H  n  C^p 

(2)  H  n  RC" 


*  <f>, 


a,  8 


(3)   At  least  two  of  the  numbers  A.,  i  =  0  ,  1 , 2 ,  .  .  .  , n, are 
of  strictly  opposite  signs. 


■1  / 


Proof:   (1)  "  (3).   The  proof  follows  from  a  sequence  of 
equivalences  below.   First,  observe  that,  by  the  last  lemma, 
H  n  Cn  g  *  ♦  is  equivalent  to  H  n  c"^*  <J>  •   Since  C^  g   is 
an  open  connected  set  in  Rn+1,  the  property  H  n  C^ &  *  <j> 
holds  if  and  only  if  H  n  co(C^g)  *  ♦.   Recall  that  co(C^) 
is  equal  to  the  open  cone  generated  by  the  row  vectors 

(qi0,qtl qin)^  i  =  °'1'2 "  °f  thG  matriX  °a,B-   HenC6' 

H  n  co(Cn  R)  *  4'  means  that  there  exist  parameters  ki  >  0, 

i  =  0,l,.'..,n,such  that  Z»=Q    C  S^i^ij  ]ijtl  '    H  °r  equivalently , 

n      n  n    n 

I    A.(  E  k.q.,)  =   E  (  E  qxiVki  =  ° 
j  =  0  3  i=0  x  1D     i-0  j=0  1D  D 


for  some  k.  >  0,  i  =  0,1,2,...  ,n.   By  lemma  2.2.1,  this  is 
true  if  and  only  if  at  least  two  of  the  numbers  Aj_  =  E-j  =  0qi:Aj' 
i  =  0,1,..., n  are  of  strictly  opposite  signs. 

(1)   «(2).   Since  C^   and  RC^   are  open  connected 
sets,  the  property  H  n  C^6  *  4>  holds  if  and  only  if 
H  n  co(C^B)  *  4-   Also,  H  n  RC^  g  *  *  if  and  only  if 
H  n  co(RC^p)  *    *.   Recall  that  co(C^)  =  cofRC^).   Hence, 
H  n  C"    M  i£  and  only  if  H  n  RCa,B  *  *' 


CHAPTER  IV 

STABILIZABILITY  OF  A  SYSTEM  WITH  DELAY 

In  this  chapter  we  wish  to  illustrate  the  applicability 
of  the  "convex  hull  technique"  which  was  used  in  the  last  two 
chapters,  to  the  study  of  the  stabilizability  of  the  following 
delayed  control  system. 

(V)       x(t)  =  ax(t)  +  bx(t-l)  +  u,  y  =  cx(t)  +  dx(t-l),   u  =  ky 

where  a,  b,  c,  d,  k  are  real  numbers. 

Our  purpose  is  to  determine  the  existence  of  a  real 
number  k  such  that  the  corresponding  characteristic  equation 
of  (P) 

A  +  (a  +  kc)  +  (b  +  kd)e~A  =  0  (4.1) 

has  only  roots  in  the  left  half  complex  plane.   Note  that 

~   ~   ~  ^2~2 

with  a,  b,  c,  d  fixed  and  c   +  d   *  0,  the  set  of  all  the 

coefficients  of  (4.1)  which  arise  from  all  possible  choices 

of  k  e  R  forms  a  straight  line 

H  =  {(a,b)|ad  -  be  -  (ad  +  be)  =  0)  (4.2) 

Now,  consider  an  exponential  polynomial  equation  of 
the  form 

A  +  a  +  be~A  =0  (4.3) 


48 


Let  Va   denote  the  set  {(a,b)  e  R  |  \    t-  a  f  be~*  =  0  implies 
ReX  <  a}.   A  hyperplane  H  is  said  to  be  V   -stabilizable  if 

H  n  Va    *    <+>.   In  particular,  when  a  =  0  we  call  0   the  domain 

of  stability  of  (4.3).   Thus,  a  system  (V)     is  stabilizable 

if  and  only  if  H  n  PQ  x    $ .   As  in  the  previous  chapters,  we 

will  study  two  main  problems  related  to  V    :   (i)  the  convex 

a 

hull  of  Pq  and  (ii)  the  characterization  of  hyperplanes 

(straiqht  lines)  which  intersect  V    . 

a 

4. 1  Convex  Hull  of  D 

2 

It  is  known  that  for  any  (a,b)  e  R" ,  equation  (4.3)  lias 

only  a  finite  number  of  roots  with  positive  real  parts 

(J.  Hale  [1],  Lemma  20.1).   This  kind  of  roots  will  be  called 

unstable  roots  in  the  following.   In  terms  of  the  D-decomposi- 

tion  method  by  Yu  Naimark  [1]  one  may  partition  the  coefficient 

2 
space  R   into  infinitely  many  regions.   We  indicate  in 

figure  4.1  these  regions  and  the  number  of  unstable  roots  of 

the  equation  when  the  coefficients  belong  to  each  region. 

The  boundary  of  these  reqions  is  the  curve  a  =  ~^  ■ ,  ■       ^ 

sin  y 

b  =  ^n  v  where  n?r  <  y  <  (n+l)7i  with  n  =  0,1,2,...  .   This 
figure  actually  suggests  most  of  the  techniques  used  in  the 
proofs  of  the  theorems  in  this  chapter. 

Hayes  [11  proved  the  following  analytic  characterization 
of  VQ. 

Theorem  4.1.1   [Hayes]  All  the  roots  of  A  +  a  +  be~A  =  0, 
where  a  and  b  are  real,  have  negative  real  parts  if  and  only  if 
(i)   a  +  1  >  0 


~b  +2  =0 


&  +  h  =0 


Figure    4 . 1 


(ii)   a  +  b  >  0 
(iii)   /a  +r   >  b   where  if  a  *  0,  r  is  the  root  of  the 
equation  r  =  -a  tan  r  in  (0,tt)  and  if  a  =  0, 
r  =  ir/2.  D 

Hence,  V n  is  the  set  of  all  points  (a,b)  which  satisfy  con- 
ditions (i) ,  (ii)  and  (iii)  of  the  Hayes  theorem.   By  the 

continuity  of  the  roots  with  respect  to  the  coefficients  of 

2 
exponential  polynomials,  one  obtains  that  PQ  is  open  in  R  . 

Also, 


Lemma  4.1.2.   P„  is  connected. 


Proof:   It  suffices  to  show  that  for  any  pair  of  points 
(alfb,)  and  (a2,b2)  form  VQ    there  is  a  path  P  c  VQ    which 
connects  them. 

Consider  the  half  line  L  =  {(a,l)|a  >  -1}.   We  assert 
that  L  c  p     Since  (a,l)  e  L  implies  a  +  1  >  0,  the  condi- 
tions (i)  and  (ii)  of  theorem  4.1.1  are  satisfied.   For 
condition  (iii),  suppose  first  that  a  *    0.   Then,  by  the 
relations  r  =  -a  tan  r  and  0  <  r  <  it  ,  one  derives 


/  2   .2  "  2 
/r  cot  r+r 


>  1 


If  a  =  0 ,  then 


77^?  =  /oTff)2  >  i  =  b 


In  both  cases  (iii)  holds.   Hence,  L  c  VQ.       Now,  for  any 
fixed  (a^b,)  and  (a2,b2),  we  construct  the  path  P  as 
P  =  P.  u  P2  u  P3,  where 


1 


(a^b)  |b  =  1  +  t  (b^l)  ,  0  <  t  <    1 


P2  =  ((.1,1)1  a  =  a1  +  t(a2-a,)l  0  <  t  <  1 
P3  =  {  (a2,b)  |b  =l+t(b2-l) ,  0  <  t  <  1} 


L  c  V, 


Note  that  if  (aQ/b0)  e  VQ,     then  for  any  b  with  -a„  <  b  <  bn, 
the  point  (aQ,b)  satisfies  (i),  (ii),  and  (iii)  of  theorem 
4.1.1.   By  this  property  one  obtains  P,  c  V  Similarly, 

P.  c  p     Hence,  P  c  p 


D 


We  remark  that  in  general  the  stability  domain  of  an 
exponential  polynomial  is  not  necessarily  connected.   For 

example,  the  stability  domain  of  x(t)  +  ax(t-l)  +  bx(t)  =  0 

2 
is  a  disconnected  open  set  in  R   bounded  by  the  straight 

n-     u         i     T\k+1  2k+l    ,  r2k+l  n2   .    n  ,  _         ,_ 
lines  b  =  (-1)     — j — Tia  +  [—  ^ — tt  ]  ,  k  =  0,l,2,...,n.   (See 

figure  4.2.)   This  suggests  an  interesting  problem  (not 

treated  here) :   "What  types  of  exponential  polynomials  have 

connected  stability  domains?" 

In  the  following  we  are  going  to  show  that 


co(PQ)  =  {(a,b)eR  |a+l>0,  a+b>0,  a-b+2>0} 


4.4) 


The  right  hand  side  will  be  denoted  by  M  in  this  chapter. 

2 
This  is  a  polygonal  convex  cone  in  R" .   To  prove  (4.4) , 

since  co(Pn)  and  M  are  open  convex  sets,  it  suffices  to 

show  that  their  closures  are  equal,  namely,  co (Fry  =  M. 

First,  we  need  the  following  lemmas. 


Lemma  4.1.3.   Let  P„  and  M  be  defined  as  above.   Then  P„  c  M. 


?-3 


Figure    4 . 2 


Proof:   Observe  that,  by  the  Hayes  theorem,  any  point  (a,b) 
chosen  from  I1   satisfies  a  +  1  >  0  and  a  +  b  >  0.   Therefore, 
one  needs  only  to  show  that  conditions  (i),  (ii)  and  (iii) 

of  Hayes  theorem  imply  a  -  b  +  2  >  0. 

7T — 2" 
Consider  a  +  2  -  /a  +r   where  r  =  a  tan  r  with  0  <  r  <  it 


when  a  *  0,  and  r=yifa=0.   Suppose  a  *  0.  Then,  by  the 
substitution  of  a  =  r  cot  r,  one  obtains 


a  +  2  -  /a2  +  r2 


/  2    2 
=  -r  cot  r  +  2  -  /r  cot  r+r 


=  lU±gos  r)  (     r    r   ,  0 
s  l  n  r         2    2 

The  last  inequality  holds,  since r- as  well  as 

tan  7-7  ace  positive  for  all  0  <  r  <  it  . 

Suppose  now  a  =  0,  in  which  case  r  =  tt/2.   Then,  one  has 


a  +  2 


rrri 

/a  +r 


=  2  - 


0 


/  2   7 

In  both  cases,  the  inequality  a  +  2  -  /a^+r   >  0  holds. 

By  (iii) ,  this  implies  that  a  +  2  -  b  >  0.   The  proof  is 
completed.  TJ 

Lemma  4.1.4.   The  point  6  =  (-1,1)  belongs  to  co  CD~^T • 

Proof:   Let  L  be  the  half  line  defined  in  lemma  4.1.2.   Recall 
that  L-,  c  P     This  lemma  follows  readily  from  the  fact  that 
for  any  c  >  0, 


(-1  +  |,1)  (    L  n  Nc(8) 


For  the  convenience   of  the  proof  of  the  next  lemma,  we 
adopt  the  following  equivalent  representation  of  M; 


M=i(a,b)  |  a  +  l>0,  b  =  ka  +  (k  +  1)   for  some  -l<k<l) 


(4.5) 


Lemma  4.1.5.   M  c  co(PQ) 


Proof:   T^et  (a0,bQ)  be  an  arbitrary  point  chosen  from  M. 
Then,  by  (4.5)  bQ  =  kQa0  +  (kQ+l)  for  some  -1  <  kQ  <  1.   We 
claim  that  for  this  fixed  kQ  there  exists  an  a  >  max{0,aQ} 
such  that  the  half  line  L   defined  as 


{ (a,b) |a>o,  b  =  kQa  +  (kQ+l) } 


lies  in  Pn .   Observe  that  if  a 


max{0,aJ,  then  for  any 


(a,b)  c    L  ,  one  has  a  +  1  >  0;  moreover, 
'       a 


a  +  b  =a  +  kQa+  (kQ  +  l)  -  (kQ  +  l)  (a+1)  >  0 

since  -1  <  k   <  1  and  a„  +  1  >  0.   Therefore,  conditions  (i) 
and  (ii)  of  theorem  4.1.1  are  satisfied  for  all  (a,b)  e    L^, 
a   being  any  number  greater  than  max{0,a0). 

For  condition  (iii) ,  it  suffices  to  show  that  for  sora 
a  >  max{0,aQ},  one  has  a2  +  r2  -  b2  >  0  f or  all  (a,b)  c  Lq, 
where  r  =  -a  tan  r  and  0  <  r  <  it.   Since  a  =  -r  cot  r  and 
0  <  r  <  tt,  it  follows  that 


2,2,2    , 
a  +r  -b   =  (-r  c 


;ot  r)2+r2-  [-kQr  cot  r  +  (kQ+D  ' 


2    2 
=  r  esc  r 

2    2  , 
=  r  esc  r I  1  - 


k»r  cot  r  -  (kg  +  1! 


(kQ  +  l)  ^J2) 


sin  r  ,  2 


Note  that  1  -  [kQcos  r  -  (kQ  +  l)    r   .1   is  continuous  in  a 
neighborhood  of  r  =  tt  ;  moreover,  one  has 


5  6 


1  -  [kQcos  tt  -  (kQ  +  i: 


sin  ii  ,  2 


1  -  k0  >  0 


Hence,  there  exists  an  e,  with  0  <  e  <  tt/2,  such  that  for  all 
r  satisfying  tt  -  z    <  r  <  n  , 
2.  2  .2 


+r  -b 


r2csc2r{l  -  [kncos  r  -  (kn  +  l)  — "^J 2  J  >°   (4-6: 


'0 


0 


Now,  let  a  =  2  max{0,a,  -  (  tt-e  )  cot  (ir-e )  }  .   Note  that  -r  cot  r 
is  an  increasing  function  of  r  in  the  interval  0  <  r  <  tt  . 
Hence,  for  every  a  >  a,  the  corresponding  r  obtained  from 
-r  cot  r  =  a  satisfies  the  inequalities  tt  -  e  <  r  <  ti .   Con- 
sequently, (4.6)  holds  for  all  (a,b)  c    LQ.   This  completes 
the  proof  of  the  claim. 

Now,  by  this  claim,  one  may  represent  (aQ,b0)  as  a  con- 
vex combination  of  the  point  (-1,1)  and  some  (a,b)  e  La  c  VQ, 
i.e.  , 


<Vb0>  -I+JT  '-1'1'  +  U  '   lTaq»(a'b)- 


Therefore,  (aQ,b0)  e  co(Pq; 


Hence,  one  concludes  that 


M  £  co(PQ) 


a-b+2>0}. 


Proof:   By  lemma  4.1.3  and  the  convexity  of  M,  one  obtains 
co(P  )  c  coTmT  =  M.   On  the  other  hand,  lemma  4.1.5  implies 


that  M  c  EoT^T.   Therefore,  M  =  co(PQ)  .   Since  for  any  m 

Y)  =  interior  (Y) 


-.m 


dimensional  convex  set  Y  c  R  ,  interior 

(see  e.g.,  Rockafellar  [ 1]  p.  46),  also  since  co(VQ)     and 


M  are  open  convex  .sets  in  R  ,  one  concludes  that.  co(f'  )  =  M  = 
{  (a,b)  |  a  +  1  >  0 ,  a  +  b  >  0,  a  -  b  +  2  >  0}  .  {} 

4 . 2  V    -Stabiliza ble  Hyper pianos 

2 

Now,  let  H  be  a  hyperplane  in  R   defined  by  Aa  +  Bb  +  C  =0. 

Since  PQ  is  an  open  connected  set,  it  follows  that  H  n  D   *  <}> 
if  and  only  if  H  n  co(PJ  *    <\> .       Observe  that,  equival  ent  ly , 
the  open  convex  cone  co ( P  )  can  be  expressed  as: 

co(PQ)  =  {(-l+k1+k2,  l+k1-k2) |k1  >  0,  k2  >  0} 

Hence,  H  n  co(P  )  *  <Ji  means  that  there  exist  k   >  0  and 
k2  >  0  such  that  (-l+k^k-,  l+k.-k  )  c  H  or  equivalently , 
A(-l+k1+k2)  +  B(l+k1-k2)  +  c  ■-=  (A+B)k.  +  (A-B)k-  +  (B-A+C)  =  0 
from  some  k,  >  0  and  k2  >  0.   By  lemma  2.2.1,  this  is  true  if 
and  only  if  at  least  two  of  the  numbers  A+B,  A-B,  B-A+C 
are  of  strictly  opposite  signs. 

Summarizing  the  analysis  above,  we  obtain  the  following 
theorem. 

2 

Theorem  4.2.1.   Suppose  H  is  a  hyperplane  in  R   defined  by 

Aa  +  Bb  +  C  =  0.  Then,  H  n  PQ  *  <Ji  if  and  only  if  at  least 
two  of  the  numbers  A-B,  A+B,  B-A+C  are  of  strictly 
opposite  signs.  [] 

Expressed  in  terms  of  the  original  parameters  a,  b,  c, 
d,  this  theorem  can  be  restated  as: 

Theorem  4.2.2.   The  system  (D)  is  V    -stabi li zabl e  if  and  only 
if  at  least  two  of  the  numbers  d-c,  d+c,  -c-d-ad-bc  are  of 
strictly  opposite  signs. 


Proof:   Recall  from  (4.2)  that  the  hyperplane  corresponding 
to  (P)  is  defined  by  da  -  cb  -  (ad+bc)  =  0.   Hence,  the 
A,  B,  C  in  the  previous  theorem  are  d,  -c  and  -(ad+bc), 
respectively.  U 

One  of  the  most  interesting  questions  in  the  stabiliza- 
tion of  a  delayed  control  system  is  the  following:   "Can   a 
system  with  delay  be  stabilized  by  a  feedback  without  delay?" 
In  general,  this  cannot  be  done  as  the  example  at  the  end  of 
this  chapter  shows.   However,  for  a  system  like  (V) ,  the 
theorem  above  tells  us  that  (V)    may  indeed  be  stabilized  by 
a  feedback  without  delay.   In  fact,  under  this  kind  of  feed- 
back, one  has  d  =  0  and  consequently,  d  -  c  =  -c  and  d  +  c  =  c 
are  of  strictly  opposite  signs. 

Given  a  e  R,  recall  that  V      is  the  set  of  all  coefficients 
(a,b)  such  that  the  real  parts  of  all  the  roots  of  X  +  a  +  be 
=  0  are  smaller  than  a.   Consider  a  substitution  of  X  by  n+ot 
in  equation  (4.3) . 

X  +  a  +  be_A  =  n  +  (a+a)  +  be~'1e~n 


The  substitution  X  =  n  +  a  induces  an  affine  transformation 

T   of  the  coefficient  space  R   of  equation  (4.3).   This 
a 

transformation  is  defined  by 

T  (a,b)  =  (x,y)  =  (a,b) 

a 


1   0 
0   c" 


+  (cx,0) 


Since  det 


1   0 


0   e 
is  defined  by 


'a  x    0,  T   is  invertible.   Its  inverse 
a 
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T"  L(x,y)  =  (a,b)  -  (x,y) 


1   0 


0   e 


+  (~a,0) 


Note  that  A  =  a  +  n  implies  that  the  condition  Re (A)  <  a  is 
equivalent  to  Re(n)<0.   Hence,  TQ   maps  Va    onto  VQ,     i.e., 

\SV  =  V 

Suppose    that    li     is    a    liyperplane    in    R      defined    by    Aa    +    Bb    + 
C    =    0 .       Then    I 
T     (II)     n    T     (V    )     -    T.  (H)     n    V, 


V      *    cf>    if    and    only    if    i,    *    T     (H      n     V  J 

u  a        a 


a  v" '       "a'     a 
Observe  that 


0 


,-1 


T  (H)  -  { (X,y) |  TuJ (x,y)  c    H} 

'1   0  )  |  .'. 

a   +  (-a,0)]W- 


=  { (x,y) |[ (x,y) 


0 


+  C  =  0} 


=  { (x,y) | Ax  +  Beay  +  C  -  Aa  =  0 } 

Therefore,  by  theorem  4.1.1  we  have  the  following  generaliza- 
tion. 


Theorem  4.2.3.   Suppose  that  H  is  the  hyperplane  in  R 
defined  by  Aa  +  Bb  +  C  =  0 .   Then  H  n  V       *  <}>  if  and  only  if 
at  least  two  of  the  numbers  A  -  Be  ,  A  +  Be    and  Be   +  C  -  A 
Aa   are  of  strictly  opposite  signs.  D 

Recall  that  for  polynomials  of  degree  n,  the  number 
inf{a|ll  n  Sn  *    4-}   is  always  bounded  from  below  for  any 
hyperplane  H.   Indeed,  theorem  2.3.1  describes  this  property. 
However  the  next  theorem  shows  that  this  is  not  the  case  for 
the  V    's. 
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2 
Theorem  4.2.4.   Suppose  H  is  the  hyperplane  in  R   defined  by 

Aa  +  Bb  +  C  =  0.   If  A  -  0,  then  inf{a|H  n  V      *    $>)    -  -». 

If  A  =  0,  then 

inf{a|HnP  *4>}  =  j nf {at R | (A-Dea) ( A+Bea) ( Bea+C-A-Aa)  =  0}  (4.7) 
1    a 

Proof:   If  A  =  0,  then  A  +  Be"   =  Bea   and  A  -  Bea   =  -Bea. 
These  numbers  are  always  of  opposite  signs  for  any  a  e  R. 
Hence,  inf{a[H  n  V      *   <j> }  =  -00.   Suppose  A  *    0.   Without  loss 
of  generality,  assume  A  >  0.   Let    denote  the  number  on  the 
right  hand  side  of  (4.7).   Note  that,  for  any  B  <  n,  the 
numbers  A  -  Be  ,  A  +  Be^   and  Be    +  C  -  A  -  A3   all  have  the 
same  sign  as  A.   Hence,  n  -  inf{a|H  n    V       *  <j>}.   On  the  other 
hand,  we  have  one  of  the  following  four  cases: 

(1)  If  n  satisfies  A  -  Be'1   =  0,  then  there  exists  an 
c  >  0  such  that  (A-BeQ)A  <  0  and  (A+Bea)A  >  0 

for  all  a    such  that  n  <  a  <  n  +  e.   Therefore 

(A-Bea) (A+Bea)  >  0  (4.8) 

for  all  a  such  that  n  <    a  <  n  +  c  • 

(2)  Suppose  that  n  satisfies  A  +  Be  '  =  0 .   Then,  as  in 
(1) ,  one  may  find  e    >    0  small  enough  such  that 
(4.8)  holds. 

(3)  Suppose  n  is  a  double  root  of  Be1  +  C  -  A  -  Aa  =  0 . 
This  implies  that  Ben  -  A  =  0.   This  case  is  then 

reduced  to  case  (1) . 

a 

(4)  If  n  is  a  root  of  multiplicity  one  of  Be   +  C  -  A  - 

Aa  =  0,  then  one  may  find  e  >  0  such  that 


r,i 


(Bt/1  +  C  -  A  -  Aa)  A  <  0  and  (A  -  Be")  A  >  0  for 
all  a  such  that  n  <  a  <  n  +  e.   Therefore, 
(Be"  +  C  -  A  -  Aa) (A  -  Be")  <  0  for  all  a  such  that 
H  <  a  <  n  +  £ .   From  the  discussion  above,  one 

concludes  that  min{a  H  n  V      *    <t> }    <  n . 

i      a 

Hence,  n  =  inf  {a|ll  n  PQ  *    4> )  .  Q 

Remark:   The  stability  domain  of  an  arbitrary   ixed  type  of 
exponential  polynomial  H(x,e  )  is  one  of  the  most  interesting 
subsets  of  the  coefficient  space  of  H(x,e' ) .   Most  of  the 
properties  of  the  stability  domains  of  polynomials  are  not 

necessarily  true  for  stability  domains  of  exponential  polyno- 

2      •,       -A 
mials.   For  example,  the  stability  domain  of  A   +  (aA  +  b) e   =0 

2 

is  the  set  bounded  by  b  =  0  and  curve  a  =  y  sin  y,   b  =  y  cos  y 

where  0  <  y  -;  -^ ■        (see  El'sgol'ts  [1])   Illustrated  by  figure 
4.3,  we  see  that  this  is  a  bounded  convex  set.   Recall, 
however,  that  the  stability  domain  of  any  polynomial  of 
degree  greater  than  or  equal  to  three  is  an  unbounded  noncon- 
vex  set.   In  chapter  one,  we  also  mentioned  that  the  stability 

domain  of  any  polynomial  is  connected.   However,  the  stability 

2       -A 
domain  of  the  equation   A   +  aAe    +  b  =  0  (see  figure  4.2) 

is  disconnected. 

It  seems  that  a  "nice"  characterization  of  the  stability 
domains  of  exponential  polynomials  is  essential  for  under- 
standing these  topological  properties  mentioned  above,  or, 
even  the  stability  domain  itself.   The  most  general  charac- 
terization available  is  a  theorem  proved  by  Pontryagin  [1] 
(see  next  section).   However,  the  difficulty  in  applying 
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T  =  o 


» 


Figure    4 . 3 
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Pontryagin's  theorem  arises  from  the  fact  that,  generally, 
it  is  not  easy  to  determine  whether  a  transcendental  equation 
has  all  the  roots  real  or  not.   Bellman  and  Cooke's  book  [1] 
contains  characterizations  of  stability  domains  for  some 
special  exponential  polynomials.   Much  still  remains  to  be 
done  for  a  better  understanding  of  the  stability  domain  of 
general  exponential  polynomials. 

4.3   An  Example 
In  this  section,  we  give  an  example  of  a  control  system 
with  delay  which  can  not  be  stabilized  by  feedbacks  without 
delay.   First,  let  us  review  a  theorem  proved  by  Pontryagin 
[1]. 


Definition:   Let  h(z,w)  =   E  a   z  w   with  m,n  nonnegative 

mn 

m,  n 

r  s 

inteqers.   A  term  a   z  w   is  called  the  principal  term  of 
^  rs 

h(z,w)  if  a    *    0   and,  if  for  each  other  term  a     with 
'        rs  mn 

a    *    0 ,  we  have  either  r  >  m,  s  >  n  or  r  =  m,  s  >  n  or 

mn 

r  >  m,  s  =  n . 


Theorem  4.3.1.   [Pontryagin]   Let  H(z)  =  h ( z , e  ),  where 
h(z,t)  is  polynomial  with  a  principal  term.   The  function 
H(iy)  is  now  separated  into  real  and  imaginary  parts  i.e., 
we  set  H(iy)  =  F(y)  +  iG(y).   If  all  the  roots  of  the 
function  H(z)  lie  to  the  left  hand  side  of  the  imaginary  axis, 
then  the  zeros  of  the  functions  F(y)  and  G(y)  are  real,  inter- 
lacing and 


G'  (y)F(y)  -  G(y)F'  (y)  >  0 


(4.9) 


(. '1 


for  each  y  c  R.   Moreover,  in  order  that  all  the  roots  of  the 
function  lie  to  the  left  of  the  imaginary  axis,  it  is  sufficient 
that  one  of  the  following  conditions  be  satisfied: 

(1)  All  the  zeros  of  the  functions  F(y)  and  G(y)  are 
real  and  interlacing  and  the  inequality  (4.9)  is 
satisfied  for  at  least  one  value  of  y; 

(2)  All  the  zeros  of  the  function  F(y)  are  real  and  for 
each  zero  y  =  yQ  of  F,  condition  (4.9)  is  satisfied 
i.e.,  F'  (y0)G(yQ)  <  0  ; 

(3)  All  the  zeros  of  the  function  G(y)  are  real  and  for 
each  zero  y  =  yQ  of  G  the  inequality  (4.9)  is 
satisfied  i.e.,  G'(y0)F(y0)  >  0. 

Lemma  4.3.2.   If  the  exponential  polynomial 

h(A)  =  eA(A2   +  PX  +  q)  +  (rX  +  m)  =  0  (4.10) 

has  only  stable  roots,  then  (q  +  m) (p  +  q  +  r)  >  0. 

Proof:   Observe  that,  by  substituting  iy  for  A,  one  obtains 
h(iy)  =  F(y)  +  iG(y)  where 

2 
F(y)  =  (q  -  y  ) cos  y  +  m  -  py  sin  y 

2 
G(y)  =  (q  -  y  ) sin  y  +  py  cos  y  +  ry 

From  Pontryagin's  [Theorem  and  the  hypothesis  that  h[X)    has  only 
stable  roots,  it  follows  that  G'(y)F(y)  -  G(y)F'(y)  >  0  for   - 
all  y  e    R.   In  particular,  if  y  =  0,  then  G'(0)F(0)  -  G(0)F'(0)= 
(q  +  m) (p  +  q  +  r)  >  0.   Now,  consider  the  following  control 
system: 


65 


E) 


X1(t) 
x2(t) 


"al  "a2 


+      u 
1 


x1(t) 

x2  (t) 


u  =  Cf  ,f  J 


xx(t) 
x0(t) 


xx(t-l) 
x2(t-l) 


where  a. ,f.;  i  =  1,2  are  real  numbers.   Note  that  there  is 
no  delay  term  in  feedback  control.   The  characteristic 
equation  of  (E)  is 

A"  +  \(a2-f2)  +  (a1-f1)  +  e~  [X-(a2-f2  +  1  +  2a1-2f1) ]  =  0 
Compare  this  equation  with  (4.10).   Then 


(q  +  m)  (q  +  p  +  r)  =  [(a^f^  -  (a2~f2  +  1  +  2a1~2f  ±)  ]  [  a,-f ,  + 


a?-f   +  1 


[ (a]-f1)  +  (a2-f2)  +  U   <  0 


for  any  real  numbers  f   and  f2>   Hence,  the  system  (E)  is 
always  unstable  for  any  feedback  without  delay. 


CHAPTER  V 
CONCLUDING  REMARKS 


Consider  a  system  (n,K)   with  paramo tors  in  K .       Tn 

chapter  2  and  chapter  3  we  discussed  respectively  the  S^- 
stabilizability  and  Cq   -stabilizabili ty  of  a  system  (it,K) 
for  which  \r.  (k)     forms  a  hyperplane.   Suppose  x^CO  contains 
only  one  point.   Then,  >;ti(K)  n  s"  *    $   is  equivalent  to 
YTi(k')  c  sn.   One  can  use  the  well  known  Routh-Hurwitz  criterion 

A  Q 

and  the  transformation  T   constructed  in  section  2.1  to  deter- 
mine if  a  point  belongs  to  Sn.   We  now  know  how  to  characterize 
S  -stabilizability  of  affine  sets  with  dimension  equal  to  n, 

Q 

n-1,  and  0.   Tor  affine  sets  with  dimensions  between  n-2  and 
1  inclusive,  no  complete  answer  is  known  yet.   Observe  that 
in  general,  for  an  affine  set  of  dimension  lower  than  n-1, 
the  fact  that  this  affine  set  intersects  the  convex  hull  of 
certain  set  S  does  not  always  imply  that  the  intersection  of 
that  affine  set  and  S  is  nonempty.   Hence,  the  "convex  hull 
technique"  gives  loss  information  for  affine  sets  with 
dimension  lower  than  n-1,  since  only  necessary  conditions  are 
obtained.   However,  because  of  their  simplicity  when  used  on 
a  computer,  those  necessary  conditions  can  be  used  as  a  quick 
check  to  eliminate  systems  which  are  "obviously"  not  Sa~stabi- 
lizable.   Skoog  [1]  obtained  some  partial  results  of  SQ-stabi- 
lizability  for  the  case  of  straight  lines  in  terms  of  Nyquist 

6  6 
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plot.   In  an  inspiring  paper,  Anderson,  Bose  and  Jury  [1] 
showed  that  Output  Feedback  Stabilization  problem  can  be 
handled  by  algorithms  from  decision  algebra  (see  Jacobson  [1]). 
This  can  be  a  very  fruitful  direction  to  pursue.   It  is  also 
interesting  to  point  out  that  the  case  of  a  hyperplane,  which 
is  the  most  complicated  among  all  affine  sets  with  dimension 
less  than  n-1  when  one   uses  decision  algebra,  turns  out  to  be 
the  easiest  for  our  approach,  but  the  case  of  a  line  which  is 
the  most  difficult  for  us  is  the  easiest  by  decision  algebra 
algori  thm. 

In  the  following,  we  will  present  a  partial  result  and  a 
computer  aided  algorithm  concerning  lines  which  intersect 
SQ.   They  should  only  be  considered  as  an  attempt  in  this 
direction . 

Proposition  5.1.1.   Suppose  f(x)  and  g(x)  are  monic  polynomials 
of  degree  n  and  n-1  respectively.   Then,  for  any  c  >  0,  A  <  0, 
there  exists  a  K  e  R  such  that  for  every  k  >  K,  n-1  of  the  roots 
of  f(x)  +  kg(x)  =  0  are  within  the  e-neighborhoods  of  the  zeros 
of  g(x)  and  the  remaining  root  is  smaller  than  A. 

Proof:   Let  A 1 ,  A2 ,  .  .  .  , Ar_ 1  be  the  n-1  roots  of  g(x)  =  0 

(counting  the  multiplicity).   For  each  i,  let  C  be  a  circle 

centered  at  Xi  with  radius  r^  <  c    such  that  the  Interior (C.) 

does  not  contain  any  A  x    \  Let  M1  =  sup  |f(z)|   and 

2         ,      ,  "    ZtCi 

Mi  "  mf  |g(z)|.   Choose  K'  >  0,  such  that  M1  <  M  K'  for  each 

zeC.  !     i 

i 

i  -  1,2,..., n-1.   Hence,  |f(z)|  <  |kg(z)|  for  any  k  >  K', 
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z  t  ci  and  i  =  l,2,...,n-l.   Therefore,  by  Rouche  theorem, 

f(x)  +  kg(x)  =  0  has  as  many  zeros  in  Interior (C . )   as  g(x)  =  0 

does  for  any  k  >  K*  and  i  =  l,2,...,n-l.   Now,  consider 

k  =  -f(x)/g(x).   Since  (-f (x) /g (x) ) '  =  -1  +  p(x)/g2(x)  with 

2 
deq  p(x)  <  deq  g  (x)  there  exist  A ±    c     R  such  that  -f(x)/g(x) 

is  strictly  decreasing  in  the  interval  (-«>,A  ).   Also,  since 

deg  f  -  deq  q  -  1,  it  follows  that  there  exists  some  A,  such 

that  k  =  -f(x)/g(x)  >  0  for  all  x  <  A2  and  lim  (-f(x)/g(x))  =  «. 

Hence,  k  =  -f(x)/g(x)  is  invertible  in  (-«,a)  with  a  =  min{A  , 
A2'A^  ~  !•   This  implies  that  for  every  k  >  -f(a)/g(a),  the  x 
obtained  from  k  =  -f(x)/g(x)  satisfies  x  <  ex  <  A ,  or  equiva- 
lently,  f(x)  +  kg(x)  has  a  zero  smaller  than  A  for  every 
k  >  -f(a)/g(a).   Now  the  proposition  follows  readily  if  one 
chooses  K  =  max{K',  -f(a)/g(a)}.  □ 

This  proposition  shows  an  interesting  geometric  property 

n  j-j 

of  SQ:   Every  straight  line  in  R   with  a  direction  (  r ,  ,  r  9  ,  .  .  .  , 

rn)  where  (r^r^  ^/r-^ r^r^  e    S*'1 ,    with  r±    *    0   inter- 

n 
sects  SQ.   Furthermore,  an  infinite  portion  of  the  line  is 

contained  in  s" 

For  straiqht  lines  which   do   not  have  such  kind  of 

"direction"  we  propose  the  following  algorithm. 

Consider  f(x)  =  xn  +  tn      a.xn~j  and  g(x)  =  Zn  nr.xn_j. 

D--L  1  -  1=1  1 

Let  f(x)  +  kg(x)  =  xn  +  l"=1(aj+krj)xn''j  =  p(x,k). 

(i)   Form  the  Hurwitz  matrix  H(k)  of  p(x,k)  treated  as 
a  polynomial  of  x. 

(ii)   Solve  the  real  roots  of  det(H(k))  =  0.   (For 
digital  computer,  this  should  not  be  a  difficult  task).   Let 
the  distinct  real  roots  be  k ,  <  kn  <  ,.  .   <    k 
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(iii)       Clioose    arbitrary    points    k.     c     (k.,k.     ,  ),   i    =    0,1,2,...,! 

1  li  +  l 

satisfying    -™    =    kn    <    kn    <    k,     <    k ,     <    k_    <    k~    <     . . .     <    k       <    k      < 
1       J  001122  mm 

km+l    =    °°* 

(iv)   Use  Lienard-Chipart  criterion  to  test  each  f  +  k.g, 

i  =  0,1,..., m.   T  f  there  exists  a  0  <  i  <  m  such  that 

f(x)  f-  k-g(x)   -  0  has  only  stable  roots,  then  for  all 

k  t  (k.,ki+1)  the  f(x)  +  kg(x)  =  0  will  have  only  stable  roots. 

Otherwise,  the  straight  line  {(a-,  +  kr,,...,a   +  kr  )  |k  e  R} 

does  not  intersect  S„. 

The  proof  of  this  algorithm  follows  from  the  next  lemma. 

Lemma  5.1.2.   Let  f(x)  and  g(x)  be  two  polynomials  of  degrees 
n  and  I    respectively  with  I    <  n.   Let  k.,  i  =  l,2,...,m  be 
the  distinct  real  roots  of  det(H(k))  =  0  where  H(k)  is  the 
Hurwitz   matrix  of  f(x)  +  kg(x).   Then,  for  each  fixed  i, 
all  the  polynomial  equations  f(x)  +  kg(x)  =  0  with  k  e     (k.,k.  ,) 
will  have  the  same  number  of  unstable  roots. 


Proof:   Let  i  be  fixed.   For  each  k  e  (k.,k.  -,),  let  p-(k), 

j  =  l,2,...,n  denote  the  roots  of  f(x)  +  kg(x)  =  0.   Recall 

that,  by  Orlando's  formula  (Gantmacher  flJ   vol.  2,  p.  196), 

det(H(k))  =  (-i)n(n+1)/2  n1?    A.(k)  *nT._(XT.  (k)  +A  (k)).   If 

j  —  i  j  r<  w   r        w 

there  exist  k'  and  k"  both  in  (k.,k.  ,-,)  and  k1  <  k",  such 

li  +  l  ' 

that  f(x)  +  k'g(x)  and  f(x)  +  k"g(x)  have  different  number 

of  unstable  roots,  then  by  the  continuity  property  of  roots 

with  respect  to  the  coefficients,  there  exists  k  with 

k.  <  k'  <  k  <  k"  <  k.+1  such  that  the  equation  f(x)  +  kg(x)  =  0 

has  either  a  pair  of  conjugate  of  purely  imaginary  roots  or  a 
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zero  root.   Hence,  det(H(k))  =  0.   This  contradicts  the  fact 

that  between  k.  and  k.,,,  there  is  no  real  root  of  det(H(k))  =0. 

Therefore,  for  every  k  c  (k.,k,+,),  the  polynomials  f  +  kg 

have  the  same  number  of  unstable  roots.  D 

By  this  1  L'inina ,  one  sees  that  it  is  enough  to  test  only 

one  point  from  each  interval  (k.  fk.  ,  ■,  )  to  determine  whether 

all  the  roots  of  f(x)  +  kg(x)  =  0  are  stable  or  not  for  any 

k  t  (k.  ,k.  ,  ,)  .   The  manor  advantaqe  of  this  algorithm  over  the 
i   l+l  J  J  J 

conventional  root  locus  method  is  that  it  involves  only  solving 
real  roots  of  an  nth  degree  polynomial.   This  is  a  well  devel- 
oped area  in  numerical  analysis.   One  may  program  that  easily 
on  a  digital  computer.   The  other  advantage  is  that  it  also 
gives  the  numerical  ranges  of  k  for  which  f  +  kg  has  only 
stable  roots.   The  major  drawback  of  this  algorithm  is  that 
it  is  very  sensitive  to  double  real  roots.   A  small  round  off 
error  may  cause  the  loss  of  a  double  real  roots  of  det(H(k))  =0. 
In  other  words,  one  may  lose  one  of  the  partition  points  k.  . 
However,  this  can  be  prevented  by  predetermining  the  number 
of  distinct  real  roots  by  Sturm's  theorem  (see  Gantmacher  CI] 
vol.  2,  p.  175).   Note  this  algorithm  also  works  for  any  poly- 
nomial of  the  form  l  .    Af.(k)x  ■*    when  f.  are  polynomials  of  k. 

3     2  2 

Example:   Suppose  f(x)  =x   +x   -x+5  and  g(x)  =  x   +  3x  -  1. 

We  wish  to  find  the  range  of  k  such  that  f  +  kg  =  0  has  only 

stable  roots. 

(i)   det(H(k))  =  3(-k  +  5)  (k  +  2)  (k-1)  . 

(ii)   The  three  roots  k=5,  k=-2,  k=l  separate  the 

real  line  into  four  intervals  (--,-2),  (-2,1),  (1,5),  (5,»). 
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(iii)   Pick  one  point  k.  from  each  interval  above  and 
substitute  them  into  H  (k)  .   It  is  easy  to  check  that  the 
only  H(k.)  which  satisfies  Lienard-Ch ipart  criterion  is  the 
one  wi  th  k .  <:  (1,5). 

Hence,  one  concludes  that  every  polynomial  equation 
f  -t  kg  --■  0  with  k  <  (1,5)  has  only  s table  roots. 

For  further  research,  this  author  feels  that  the  study 

of  the  S  -stabi lizabi 1 ity  of  affine  sets  with  dimension 
a 

strictly  lower  than  n-1  deserves  an  immediate  attention. 

Besides  the  engineering  application,  this  study  also  provides 

a  better  understanding  of  the  set  S   itself.   After  this,  the 

a 

investigation  of  S  -stabi lizability  of  systems  with  multiple 
control  would  be  the  next  goal.   For  nonlinear  control  systems, 
a  problem,  initially  proposed  by  Letov  ill,  of  determining  the 
minimal  number  state  measurements  needed  for  stabilizing  the 
system  by  a  feedback  control  was  considered  by  Casti  and 
Letov  [ 1] . 

In  a  practical  situation,  it  is  not  enough  just  to  know 
that  the  system  is  stabilizable  by  a  constant  gain  feedback. 
The  numerical  scheme  for  finding   a   feedback  gain  matrix  is 
also  important.   The  papers  by  Miller,  Cochran  and  Howze  CI]; 
Luus  Til,  McBrinn  and  Roy  Cl],  as  well  as  Sirisena  and  Choi  [1] 
represent  some  of  the  research  in  this  direction. 


APPENDIX    I 


G  i  veil    ex    i     l<    ruca  J  1    t  ha  t 


S       =    I  (  a  ,  ,  a .,  , 
a  12 


.  ,  a    )     t     R         x       +    L.     ,a.x       J    -    0    implies 


Re(x)     <    a; 

In  this  appendix  we  wish  to  show  that  S   is  connected.   Since, 
in  Rn ,  connectedness  is  equivalent  to  path  connectedness,  it 

,  a 


suffices  to  show  that  for  every  a'  =  (aJ,ai, 


e    Sn, 


a"  =  (aV  ,  a  ~  ,  .  .  .  ,  a"  )  e  S   there  exists  a  path  p:  [0,1]  ■*  R 
1  '  2       n      a 

such  that  p(!0,ll)  <=  Sn  with  p(0)  =  a'  and  p(l)  =  a".  Now 
consider  a  mapping  ti  from  C  into  C  where  C  is  the  set  of 
all  complex  numbers,  such  that  for  each  (x, , x0 , . . . , x  )  e  C  , 


i!  (x,  ,x2,  .  .  .  ,  Xj 


(al'a2 an} 


where  a.'s  satisfy  the  relation  {x!  x   +  Z.  ,a.x   J  =  0}  = 

(x,  ,xn,  ...  ,x  1.   Note  that  ir  is  continuous.  Also,  if 
1  '2.                n 

(x,,x9,...,x  }  is  symmetric  with  respect  to  the  real  axis 

i.e.,  {x1,x~,....x  }  =  {x,,x.,,...,x  }  where  x-  represents  the 

'    1 '  2       n       1  '  2 '    '  n  j    l 

coniuqate  of  >:.,  then  n  (x-,  ,  x„  ,  .  .  .  ,  x  )  e  R  . 
j  j         j  1 '  2       n 

Now,  for  the  a'  and  a"  given  above,  arrange  the  roots  of 


n     n    ,  n-j 
x   +  l  .     ,  a  '.  x 

1  =  1  1 


,n 


=  0andx   +  Z  .     ,a'.'x 
1  =  1  ] 


-  n~l  - 


=  0  such  that 
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<*|     ^=1a^n"j)    =    (x-    x- *Vx2m+l >:;} 


U|    I 


I"1  J 


,11           „     ii-  J  i            r  vn     vn  vn        v  ■• 

'j=lajx         }    =    ^xrx2 x2r'X2r 


+  1'  * 


where  x!,  i  =  2m+l,...,n  and  x",  j  =  2r+l,...,n  are  real  roots; 


x '  ,  1 

l 

such  that  x 


12  ...,2m,  and  x '.  ,  j  =   l,2...,2r  are  complex  roots 


A2i-1 


,  3 


=  1  ,  2 ,  .  .  .  , m  and  x " 


n 


x3j-ll  j  =  1,2, 


,r.   Note  that  i(x|,xJ,...,xM  =  (aj_,  .  .  .  ,a^)  and  tt(x'^,x£, 


•  'x. 


.  a ,  ,  a  0  , 


n-    '  i'":. 
defined  by 


,a").   Consider  the  path  p:   [0,1]  +  C 


x')  +  (1-t) (x",x"  . . . ,x") 

n  1   ^       n 


p(t)  =  t(x^,x^, 

Note   for  each  t  c    [0,1]  the  set  {tx[  +  (l-t)x^,  txj  +  (l-t)xj, 
...,tx'  +  (l-t)x"},  is  symmetric  with  respect  to  the  real 


n 


.  n 


axis.   Hence,  tt  (p  ( !.  0  ,  1  ] )  )  -    R"  .   Also,  since  a'  c    S^    and 
a"  £  sn,  by  the  definition  of  sJJ,  one  has  Re(x!)  <  a   and 


a 


Re(xn )  <  a   for  every  i  =  1,2, ...,n.   It  follows  that  for  any 
t  e  [0,1]  and  for  every  i  =  1,2, ...,n, 

Re(tx!_  +  (l-t)xV)  =  t  Re(x[)  +  (l-t)Re(x?)  <  a 

Therefore,  ir  (p  ([  0, 1])  )  c  s".   Note  that  ir(p(0))  =  a"  and 
tt(p(D)   =  a'.   Hence,  tt  (p)  is  a  path  in  Sn  which  connects 
a'  and  a".   This  proves  s"  is  a  connected  set  in  R  . 
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For  a  lixcil  n  let  q"  (,  »  (q^)  bo  Llic  matrix  roprcauntn- 

tion  of  T     defined  in  section  3.1.   Recall  that  for  each 
a ,  (5 


fixed  i,  the  row  (q^^^ qin)  of  Qa ,  g  reP 

coefficients  of 

n 


resents  the 


:s-a)n-1(s-B) 


=  0 


q  .  .  s 


(II. 1) 


Now,  let  f(x,y)  =  x"  1y1.   Then, 


/       x  "_1  /     nil 

f(-a+s,-B  +  s)  -=  (s-a)     (s-h) 


11.21 


In  terms  of  Taylor's  theorem,  one  has 
f(a+h,  b  +  k)  =  (n-j)  ! 


1     Z  (h*  +  k^)n^(xn-V) 


j  =  0 


'  9  x     d  y 


x=a 

y  =  b 


By  substituting  a  =  -a,  b  =  -B  and  h  =  k  -  s  into  the  above 
equation  one  obtains 


n-i  ,   «»i 


By  comparing  the  right  hand  sides  of  (II. 1)  and  (II. 3),  one 
derives  another  expression  for  the  entries  of  QQ  „,  namely, 


n_:>   (II. 3) 


'ID 


1     <±    +  A)n-3(xn--  ' 


-j)  !  l3x    3y 


y 


x=-a 
y=-P 


(II. 4) 
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n       ,   o  v  n  (n+1)  /2 

To  show  that  for  each  fixed  n,  det(Qa>e)  -  (a-fcU 

let  Mn(x,y)  =  (^n  •)  be  the  (n+1)  x  (n+1)  matrix  with 


m".  o  _^-  (4+  ^)n"3(xn-V) 


ij  "'  (n-j)  !   3x    dy 
We  claim  that  det(Mn(x,y))  -  (y-x)  "  (,,  +  1) /2  .   To  prove  this, 
we  use  mathematical  induction  on  n.   Let  M  =  (mjfj)  be  the 
matrix  obtained  from  Mn(x,y)  by  starting  from  the  (n+l)-th 
row  of  Mn(x,y),  subtracting  successively  the  (j-l)-th  row 

from  the  j-th  row,  j  =  n  +  l,n 2.   These  are  elementary 

row  operations  which  do  not  change  the  determinant  of  Mn(x,y). 
Therefore,  det (Mn(x,y) )  =  det(M).  Furthermore,  for  any  i  >  2 
3  ^  2, 

A)n-J[(xn-1yi"1)(y-x)J 


(n-j) !  3x    dy 


+  a^u^y-1)  +  (^"V"1) 


(3  +  3)n^(y-x)} 

d  x    a  y 


JL^-i   Iv"-1,,1"1' 
(n-j).l  ^  *'  v3x  '  dy' 


(y_x)(JL  +  ^)"-J(x"  V 


.    .  n-1 

(y-x)mi-i,j-i 


Also,  m.,  =  0  for  i  >  2,  since  m^  =  1  for  each  i  =  1,2,..., 
n+1.   Hence,  by  factoring  (y-x)  out  of  2nd,  3rd,  and  (n+l)-th 
row  of  M  and  applying  induction  hypothesis,  one  obtains 
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det(Mn(x,y))  =  (y-x)ndet(Mn~1(x,y))  =  (y-x)  n  •  (y-x)  (n~  l)  n/2 
=  (y-x)n(n  +  1>/2 

How,  observe  that  Mn(-a.-B)  =  Qn  „ .       Therefore 

a,  3  ' 

-a ,  p.  ■ 


detect   )  =  dct(Mn(x,y)) I  _    =  (a-B)n(n+1)/2 


y=-B 
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